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SECTION  1 
INTRODUCTION 


•i. 

The  HFNET  computer  code  is  a  simulation  tool  for  estimating  nuclear 
effects  on  HF  communications  systems.  The  code  has  been  developed  under 

v 

Defense  Nuclear  Agency  (DNA)  sponsorship  for  the  purpose  of  rapid  analysis 
of  nuclear  effects  on  multiple  HF  to  VHF  radio  links  in  a  multiple  nuclear 
burst  environment .  The  basic  philosophy  behind  HFNET  has  been  to  include 
all  effects  believed  to  be  important  to  HF/VHF  communications,  but  to  do  so 
in  a  manner  which  gives  adequate  (not  necessarily  the  best)  accuracy  while 
avoiding  long  running  times  or  large  memory  requirements. 

In  most  cases  the  goal  of  producing  a  fast  running  code  has  re¬ 
sulted  in  simplified  models.  For  example,  the  nuclear  environment  models 
in  HFNET  are  much  less  elaborate  than  those  found  in  detailed,  state-of-the- 
art  numerical  simulations  such  as  MICE  or  MELT  (Reference  1) .  The  HFNET 
phenomenology  models  use  orders  of  magnitude  less  computer  resources  while 
still  providing  a  time -dependent  specification  of  nuclear  debris  parameters 
which  is  adequate  for  HF  propagation  effects  calculations. 

On  the  other  hand,  the  goal  of  being  able  to  handle  large  numbers 
of  links  and  bursts  has  forced  a  fairly  elaborate  simulation  structure  and 
sophisticated  data  management  techniques  to  be  incorporated  in  the  code. 

These  features,  along  with  the  fast  running  models,  enable  HFNET  to  handle 
problems  much  too  large  for  other  nuclear  environment  HF  communications  codes ^ 

This  report  describes  the  September  1979  version  of  HFNET  and  is 
primarily  oriented  towards  the  user  of  the  program.  As  such,  the  most 
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detailed  information  concerns  the  input  and  output.  No  attempt  has  been 
made  here  to  fully  explain  the  various  physical  and  mathematical  models 
incorporated  in  the  code.  Any  discussion  of  these  models  is  made  in  the 
context  of  explaining  some  facet  of  the  input  or  output.  Full  descriptions 
of  the  models  in  HFNET  are  contained  in  a  companion  report  (Reference  2) . 

Sections  2  and  3  to  follow  describe  the  input  requirements  and  out¬ 
put  options  available.  The  input  section  gives  details  on  the  format,  units, 
defaults  and  restrictions  of  the  various  input  quantities  as  well  as  examples 
of  their  use.  The  output  section  "walks  through"  the  output  from  an  elaborate 
sample  problem  and  explains  in  detail  the  quantities  displayed  in  the  various 
output  formats. 

The  information  contained  in  the  three  appendices  is  oriented 
more  towards  the  implementor  of  the  code  than  the  general  user.  Appendix 
A  gives  information  on  the  simulation  structure  used  by  HFNET.  Anyone  not 
familiar  with  event -structured  simulation  techniques  should  read  this  appen¬ 
dix  to  get  a  feel  for  how  the  program  works.  Appendix  B  contains  alphabet¬ 
ical  lists  of  the  routines  and  common  blocks  in  HFNET  with  short  summaries 
of  how  each  is  utilized  by  the  code.  Finally,  Appendix  C  contains  a  short 
discussion  of  the  implementation  of  the  code  from  a  computer  programmer's 
perspective. 
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SECTION  2 
HFNET  INPUT 


2.1  GENERAL 

Input  to  HFNET  comes  from  two  major  sources:  a  user -supplied  input 
file  containing  problem  specification  and  program  control  information,  and 
several  binary  data  files  containing  fixed  databases  required  by  the  program. 
The  user-supplied  input  file  is  the  subject  of  this  section  while  the  binary 
data  files  are  discussed  in  Appendix  C. 

The  input  file  consists  of  a  sequence  of  standard  (up  to  80  charac¬ 
ter!  input  records,  or  lines.  These  input  lines  are  organized  into  groups, 
depending  on  the  type  of  information  they  contain,  with  each  group  identified 
by  a  special  alphanumeric  identifier  called  a  "keyword".  Following  the  key¬ 
word  line  are  one  or  more  data  lines  whose  content  and  format  vary  from  key¬ 
word  to  keyword.  In  addition  to  keyword  and  data  lines  the  input  file  may 
also  contain  any  number  of  comment  lines. 

HFNET  problem  specification  and  simulation  control  information  are 
read  at  the  beginning  of  each  run  by  the  input  module.  The  input  module  does 
more  than  "just  read  input",  it  translates  the  external  problem  specification, 
contained  in  the  input  file,  into  an  internal  specification,  ready  for  execu¬ 
tion  by  the  simulation.  During  this  translation  process  the  input  module 
performs  several  tasks: 

1)  reads  input  lines  and  saves  them  for  later  echoing  on  the 
output , 

2)  recognizes  keyword  lines  and  checks  for  keyword  errors, 

3)  interprets  data  lines  depending  on  the  keyword,  supplies 
default  values  when  needed,  and  checks  for  data  errors. 


4)  skips  comment  lines. 


5)  stores  events  in  the  event  list, 

6)  stores  data  in  dynamic  storage  and  common  blocks,  and 

7)  causes  the  program  to  abort  if  input  errors  are  detected. 

As  can  be  seen  above,  many  checks  are  made  by  the  input  module  to  insure 

that  the  input  is  "reasonable".  This  results  in  many  runs  being  aborted 
before  they  even  get  started,  usually  due  to  some  trivial  error  in  the  in¬ 
put.  The  desirability  of  this  program  "feature"  is,  hopefully,  self  evident. 

Keyword  lines  are  used  to  inform  the  input  module  what  kind  of 
data  line(s)  follow.  A  keyword  line  contains  one  of  the  input  keywords 

listed  in  Table  1.  A  few  simple  rules  govern  the  use  of  keywords: 

1)  keywords  must  be  left  justified  on  the  keyword  line, 

2)  except  for  $RUN,  no  keyword  is  required, 

3)  except  for  $RUN,  keywords  may  appear  in  any  order, 

4)  the  $RUN  keyword  is  required  and  must  be  last, 

5)  no  keyword  may  appear  more  than  once. 

If  any  of  these  rules  are  violated  the  run  will  be  aborted  with  an  error 
message  explaining  the  problem. 

Following  each  keyword  are  one  or  more  data  lines.  Data  lines 
may  contain  either  numeric  or  alphanumeric  data.  Data  lines  containing 
numeric  data  are  read  using  FORTRAN  list-directed  ("*")  format  into  float¬ 
ing-point  (REAL)  variables.  In  practice  what  this  means  is  that  numeric 
data  items  may  be  expressed  in  any  "natural"  format  (with  or  without  deci¬ 
mal  points  or  E's)  and  they  may  appear  anywhere  on  the  input  line,  separated 
by  any  number  of  spaces  or  tabs  and/or  a  single  comma.  Alphanumeric  data 
lines  are  read  using  FORTRAN  "A"  format  and  must  be  left  justified  on  the 
input  line.  The  quantities  contained  on  data  lines  and  the  order  in  which 
data  lines  must  appear  are  explained  below  in  the  subsections  devoted  to 
the  individual  keywords. 
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Table  1.  HFNET  keywords. 


Keyword 

Description 

$ I DENT 

Problem  identification  lines 

$AMBIENT 

Ambient  problem  conditions 

$BURSTS 

Nuclear  bursts 

$CLD0UT 

Cloud  list  output 

$N0DES 

Define  nodes  (transmitters/receivers) 

SLINKS 

Define  links 

$HFCALC 

HF  propagation  calculations 

$PLCALC 

Plume  mode  propagation  calculations 

$DEBUG 

Debug  output 

$RUN 

Simulation  start/stop 
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Any  input  line  which  is  either  totally  blank  or  begins  with  a  semi¬ 
colon  is  treated  as  a  comment  line.  Comment  lines  are  echoed  on  the  output 
but  are  otherwise  ignored  by  the  input  module.  There  are  no  restrictions  on 
the  use  of  comment  lines;  any  number  of  comment  lines  may  appear  anywhere  in 
the  input  file.  The  user  would  be  well  advised  to  take  full  advantage  of  the 
power  of  eminent  lines  to  document  exactly  what  problem  is  being  run  (and, 
perhaps,  wny  the  problem  is  being  run).  Comment  lines  are  also  extremely  use¬ 
ful  as  labels  for  input  quantities  as  shown  in  the  examples  given  below. 

HFNET  keeps  track  of  time,  simulation  time,  in  units  of  "seconds 

from  midnight  (GMT)  on  day  zero".  Day  zero  is  the  date  (year,  month,  day) 

specified  on  the  $AMBIENT  data  line.  In  these  units,  noon  on  day  zero  is 

t  .  =  43,200  seconds  and  19:00  on  day  three  is  t  .  =  327,600  seconds. 

Needless  to  say,  this  method  of  timekeeping  is  fine  for  computers  but  not 

for  people.  In  an  attempt  to  make  time  easier  to  express,  HFNET  allows  the 

user  to  input  time  in  a  special  seven  digit  format,  "DHHMM.SS",  called  DBMS 

format.  In  this  format  noon  on  day  zero  is  t_.„w„  =  01200.00  =  1200,  and 

19:00  on  day  3  is  =  31900.00  =  31900.  Both  absolute  (GMT)  time  and 

time  intervals  are  expressed  in  DHMS  format.  For  example,  a  time  interval 

of  45  minutes  is  dt  =  00045.00  =  45.  The  components  (D,  HH,  MM  and  SS)  of 

DHMS  format  are  restricted  as  follows:  .1  —  Ds9,  0;SHH^24,  0^MM^60  and  O^SS^GO. 

If  these  restrictions  are  violated  the  program  will  abort.  An  example  of  bad 

DHMS  format  is  dtnmir  =  90  (used,  presumably,  to  indicate  90  minutes) .  The 
UHM5 

"correct"  way  to  specify  90  minutes  is  dt^^g  =  130. 

The  following  ten  subsections  contain  discussions  of  the  individual 
keywords,  their  use  and  their  data  requirements.  In  addition,  Tables  2 
through  11  give  a  "quick  look"  summary  of  each  keyword  including  examples 
of  their  use.  These  tables  contain  a  few  notational  oddities  used  to  indi¬ 
cate  the  range  of  acceptable  values  a  given  data  item  may  take  on.  The 
notation  [a,b]  is  used  to  indicate  that  the  associated  data  item,  say  X, 
must  satisfy  the  inequality:  a^X^b.  A  "U"  is  used  to  denote  the  set 
theoretic  union  of  two  or  more  such  intervals.  For  example,  [0]U[2,30] 


should  be  translated  as  X  =  0  OR  2  £  X^30  ,  and  [0,+]  translates  into 
X  ^0  .  The  input  module  checks  all  input  data  against  the  specified  limits 
and  aborts  the  program  if  any  data  item  is  found  to  be  out  of  range. 

2.2  $IDENT  KEYWORD 

The  $IDENT  keyword  is  used  to  input  arbitrary  problem  identifica¬ 
tion  information  that  will  be  echoed  on  the  summary  output  as  a  means  of 
easily  identifying  the  problem  being  run.  Any  number  of  problem  identifica¬ 
tion  lines  may  appear,  each  containing  up  to  80  arbitrary  characters. 

2.3  $AMBIENT  KEYWORD 

The  ambient  "weather"  conditions  for  a  problem  are  specified  using 
the  $AMBIENT  keyword.  The  quantities  associated  with  the  $AMBIENT  keyword 
are  the  date  (year,  month,  day),  the  sunspot  number,  the  magnetic  storm  in¬ 
dex,  the  solar  disturbance  index,  the  east  and  north  wind  velocities  and 
a  wind  model  selector.  These  quantities  are  all  expressed  on  a  single  data 
line  which  follows  the  $AMBIENT  keyword  line.  If  $AMBIENT  data  are  not 
specified  the  defaults  shown  at  the  bottom  of  Table  3  are  used. 

The  values  specified  on  the  $AMBIENT  data  line  are  used  by  the 
different  ambient  environment  models  in  various  ways.  The  current  version 
of  the  code  does  not  use  the  year  parameter  at  all.  The  month,  day  and 
sunspot  number  are  used  by  all  three  ionospheric  models  and  by  other  models 
such  as  the  ambient  D-region  model.  The  magnetic  storm  index  is  used  by 
the  RADC  (polar)  ionospheric  model  exclusively.  Likewise,  the  ambient  D- 
region  absorption  model  is  the  only  model  which  uses  the  solar  disturbance 
index.  Finally,  the  three  wind  parameters  are  used  to  select  one  of  three 
wind  models  (constant  winds,  median  winds  or  diurnal  winds)  which  affect 
the  behavior  of  late  time  nuclear  debris  clouds. 
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Table  2.  $IDENT  keyword. 


Keyword  line:  $ I DENT 

Data  line:  HIDENT 


Data  name  Description 

HIDENT  Problem  identification  line  (80  characters) 


Example 


:  $ I DENT  EXAMPLE 

;  PROBLEM  IDENTIFICATION  LINES  FOR  AN  EXAMPLE  HFNET  PROBLEM. 
SI DENT 


HIDENT 


EXAMPLE  HFNET  PROBLEM  FOR  THE  USER'S  MANUAL 


THF  PROBLEM  CONSISTS  OF  3  BURSTS.  4  NODES.  3  LINKS,  CALCULATIONS  EVERY  1/2 
HOUR  FOR  2  HOURS  AND  EVERY  HOUR  THEREAFTER  FOR  1  DAY.  STARTING  AT  19:00  (GMT) . 


Notes: 


(1)  Any  number  of  $IDENT  data  lines  may  appear. 

(2)  Problem  identification  lines  are  echoed  on  the  first  page 
of  the  summary  output. 
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Table  3.  $AMBIENT  keyword. 


Keyword  line:  $AMBIENT 

Data  line:  YEAR,  MONTH,  DAY,  SSN,  KP,  SOLAR,  VEWIND,  VNWIND,  WINDS 


Data  name 


Description 


YEAR 

MONTH 

DAY 

SSN 

KP 

SOLAR 


VEWIND 

VNWIND 

WINDS 


Year  [1900,  2000] 

Month  [1,12] 

Day  of  the  month  [1,  31] 

Zurich  smoothed  sunspot  number  [10,  200] 
Magnetic  storm  index  [0,  9] 

Solar  disturbance  coefficient  [0,  2] 

0  =  quiet 

1  =  disturbed 

2  =  highly  disturbed 

East  wind  velocity,  wind  model  0  (m/s) 

[-100,  100] 

North  wind  velocity,  wind  model  0  (m/s) 

[-100,  100] 

Wind  model  selector  [0,  2] 

0  =  constant  winds  (VEWIND/VNWIND) 

1  =  median  winds 

2  =  full  diurnal  winds 


Example 


;  $AMBIENT  EXAMPLE 

;  DATE  =  1  APRIL  1984  —  SUNSPOT  NUMBER  =123  --  KP  INDEX  =  4 

;  SOLAR  =  QUIET  —  WIND  MODEL  =  MEDIAN  bJINDS 


$A MB  I ENT 

■•'.’EAR  MONTH  DAY  SSN  KP  SOLAR  VEUINP  VNU1NB  LUNDS 

1984  4  1  113  4  0  0  0  1 

Notes:  (1)  At  most  one  $AMBIENT  data  line  may  appear. 

(2)  Defaults  for  $AMBIENT  data  are:  YEAR  =  1979;  MONTH  =  1; 
DAY  >  1;  SSN  -  100;  KP  *  2;  SOLAR  =  VEWIND  =  VNWIND 
=  WINDS  =  0.  These  values  are  used  when  no  $AMBIENT 
data  are  supplied  in  the  input  file. 
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2.4 


$BURSTS  KEYWORD 


Nuclear  bursts  are  specified  using  the  $B(JRSTS  keyword.  Any 
number  of  bursts  may  be  specified  with  almost  any  combination  of  the  follow¬ 
ing  burst  parameters:  time,  latitude,  longitude,  altitude,  yield,  fission 
yield.  X-ray  yield,  late-time  radial  spread  rate,  conjugate  burst  flag 
and  output  flag.  However,  due  partly  to  uncertainties  in  phenomenology  and 
partly  to  modeling  inadequacies  there  are  a  few  arbitrary  restrictions  on 
the  location  and  yield  of  bursts.  These  restrictions  include: 

•  the  location  of  a  burst  (lat.  Ion)  cannot  be  nearer  than 
one  degree  from  any  pole  (geographic  or  geomagnetic) , 

•  the  burst  altitude  must  be  in  the  range  0  to  1000  kilometers 
(no  underground,  underwater  or  outer  space  bursts  allowed), 

•  the  total  yield  must  be  at  least  1  kt  (.001  NTT)  and  no  more 

than  10*-2"  This  odd  restriction  translates  into 

maximum  yields  of  100,  56,  10,  1  and  .32  megatons  at  alti¬ 
tudes  of  0,  100,  400,  800  and  1000  kilometers  (respectively) 

As  usual,  the  program  will  abort  if  any  burst  is  specified  with  parameters 
which  violate  these  restrictions. 

Late  time  wind  dispersal  of  subsiding  nuclear  debris  clouds  can 
be  simulated  in  two  ways  within  HFNET.  The  most  straightforward  way  is  to 
select  one  of  the  internal  wind  models  (using  the  $AMBIENT  keyword)  and 
let  the  wind  blow  the  debris  where  it  may.  The  other  way  is  to  use  the 
late  time  radial  spread  rate  parameter  on  the  $BURSTS  data  line  to  force 
the  debris  to  expand  at  a  constant  rate  with  no  associated  movement  of  the 
center  of  mass.  The  two  approaches  can  also  be  combined,  with  the  effects 
simply  added  together. 

To  more  realistically  model  a  very  high  altitude  explosion,  the 
debris  and  energy  is  split  into  two  components.  One  component  moves  down 
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Table  4.  $BURSTS  keyword. 


Keyword  line: 
Data  line: 


$BURSTS 

TIME,  LAT,  LON,  ALT,  Y,  FY,  XY,  SPREAD,  CBFLAG,  PBFLAG 

Data  name  Description 

TIME  Burst  time  ( DHMS )  [0,  +] 

LAT  Burst  latitude  (®N)  [-89,  89] 

LON  Burst  longitude  (°E)  [-360,  360] 

ALT  Burst  altitude  (km)  [0,  1000] 

Y  Yield  (MT)  [.001,  10(2"  Alt/400)] 

FY  Fission  yield  (MT)  [0,  Y] 

0  =  use  FY  =  (Y  if  Y  < .01 
] Y/2  if  Y^  .01 

XY  X-ray  yield  (MT)  [-1]  U  [0,  Y] 

0  =  use  XY  =  3/4  Y 

-1  =  use  XY  =  0 

SPREAD  Late  time  radial  spread  rate  (m/s)  [0,  1000] 

CBFLAG  Conjugate  burst  flag  [0,  1] 

0  =  suppress  conjugate  bursts 
1  =  allow  conjugate  bursts 

PBFLAG  Print/binary  output  flag  (packed,  PB)  [00,  21] 

P:  0  =  off 

1  =  on,  output  at  cloud  creation  and 

spl its 

2  =  on,  output  at  cloud  creation,  splits 

and  updates 

B:  0  =  off 

1  =  on,  binary  output  at  cloud  creation 
and  updates 
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Table  4.  (cont.) 


Example 


;  $0URSTS  EXAMPLE 


;  SETUP 

THREE 

BURSTS: 

(1)  1 
(2)  20 
(3)  5 

MT,  200 
MT.  1 
MT,  10 

KM.  40 
KM,  40 
KM,  40 

DEG 

DEG 

DEG 

N.  90 
N,  110 
N,  70 

DEG  U. 
DEG  LI, 
DEG  U, 

18:06  GMT 
18: 14  GMT 
18:31  GMT 

SBURSTS 

:  T I  ME 

LAT 

LON 

ALT 

Y 

FY 

XY 

SPREAD 

CBFLAG 

PBFLAG 

01806 

40 

-90 

200 

1 

.  5 

.75 

0 

0 

10 

01814 

40 

-110 

1 

20 

10 

15 

0 

0 

00 

01831 

40 

-70 

10 

5 

2.5 

3.75 

0 

0 

00 

Notes:  (1)  Any  number  of  $BURSTS  data  lines  may  appear 

(2)  Notice  the  restrictions  on  burst  location  and  yield; 

the  yield/altitude  restrictions  result  in  the  following  table: 


Altitude  (km) 
0 

100 

400 

800 

1000 


Max,  yield  (MT) 
100 
56 
10 
1 

.32 


(3)  The  output  flag  only  affects  the  single  burst  being  set  up; 
output  for  the  entire  cloud  list  can  be  independently  speci¬ 
fied  using  the  $CLD0UT  keyword. 
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the  local  magnetic  field  line  until  it  encounters  enough  air  to  contain  it, 
then  acts  as  if  an  explosion  with  appropriately  reduced  yield  occurred  at 
the  displaced  position.  The  other,  usually  much  smaller,  component  moves 
to  the  magnetic  conjugate  region  and  acts  as  if  it  were  an  additional  ex¬ 
plosion  there.  The  conjugate  (pseudo)  explosion  occurs  with  a  delay  time 

g 

appropriate  to  a  debris  velocity  of  10  cm/sec  to  the  conjugate  region. 

The  conjugate  burst  flag  allows  the  user  to  suppress  the  formation  of  this 
conjugate  burst.  Turning  off  conjugate  bursts  is  advisable  when  the  effects 
of  debris  in  the  magnetic  conjugate  region  would  be  negligable. 

The  print/binary  output  flag  controls  the  outputting  of  cloud 
parameters  to  the  printed  output  file  and  to  an  auxiliary  binary  data  file 
used  to  generate  plots.  Three  levels  of  printed  output  are  available: 

(1)  none,  (2)  at  cloud  creation  and  split  times,  and  (3)  at  cloud  creation 
split  and  update  times.  The  user  should  be  warned  that  the  amount  of  prin¬ 
ted  output  generated  at  cloud  update  times  (level  3  output)  can  be  excessive. 
The  binary  (plot)  output  is  primarily  intended  for  debugging  the  nuclear 
phenomenology  models  and  requires  an  additional  program  to  process  the  data 
into  plots. 

2.5  $CLD0UT  KEYWORD 

Output  from  the  nuclear  phenomenology  models  can  be  requested  in 
two  different  ways.  As  explained  in  the  previous  section  the  output  flag  on 
the  $ BURSTS  data  line  can  be  used  to  request  output  describing  individual 
clouds  at  cloud  creation,  split  and  update  event  times.  Independently  the 
$CLDOUT  keyword  can  be  used  to  request  that  data  describing  every  cloud  in 
the  cloud  list  be  output  at  simulation  times  specified  by  the  user.  In  most 
cases  the  $CLDOUT  keyword  method  is  preferable  since  it  can  provide  "snap¬ 
shots"  of  the  nuclear  environment  at  precisely  the  same  times  that  HF  propa¬ 
gation  and/or  plume  mode  calculations  are  made.  This  can  be  very  useful 
when  the  results  of  these  calculations  are  later  analyzed. 
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Table  5.  $CLDOUT  keyword. 


Keyword  line:  $CLD0UT 

Data  line:  Tl,  T2,  DT,  PBFLAG 


Data  name 

')  ^.ription 

\ 

Tl 

A 

First  event  time  (DHMS)  [0,  +] 

T2 

Last  event  time  (DHMS)  [0,  +] 

DT 

Event  rescheduling  interval  (DHMS)  [0,  +] 

PBFLAG 

Print/binary  flag  (packed,  PB)  [00,  11] 

P:  0  =  off 

1  =  on,  printed  cloud  list  output 

B:  0  =  off 

1  =  on,  binary  cloud  list  output 


Example 


SCLDOUT  EXAMPLE 


;  CLOUD  OUTPUT  EVERY  HALF  HOUR  FROM  19:00  TO  21:00  (GMT) 
;  AND  EVERY  HOUR  FPOM  21:00  TO  19:00  (THE  NEXT  DAY). 


SCLDOUT 


a  Tl 

T2 

DT 

PBFLAG 

01900 

02100 

0030 

10 

02200 

11900 

0100 

10 

Notes:  (1)  Any  number  of  $CLD0UT  data  lines  may  appear. 

(2)  Cloud  list  binary  output  goes  to  an  auxiliary  output  file 


The  $CLDOUT  keyword  is  used  to  set  up  cloud  list  output  events. 

The  $CLDOUT  keyword  line  can  be  followed  by  any  number  of  data  lines  each 
containing  four  quantities:  first  event  time,  last  event  time,  event  resched¬ 
uling  interval,  and  print/binary  output  flag.  The  printed  and  binary  out¬ 
put  are  independently  controllable  using  the  output  flag  and  each  contains 
detailed  information  describing  every  cloud  in  the  cloud  list  at  the  time 
of  the  event.  The  printed  output  is  written  to  the  standard  output  file 
and  the  binary  output  to  an  auxiliary  output  file. 

In  addition  to  the  detailed  printed  and  binary  output  described 
above,  cloud  list  summary  output  is  also  generated  at  each  cloud  list  out¬ 
put  time.  This  output  is  not  controlled  by  the  print/binary  output  flag. 
Examples  of  cloud  list  output  are  given  in  Section  3. 

2.6  $N0DES  KEYWORD 

In  the  HFNET  code  a  node  can  be  a  transmitter  or  a  receiver  or 
both.  Nodes  are  created  using  the  $NODES  keyword  and  then  connected  into 
links  (transmitter/ionosphere/receiver  triplets)  using  the  $LINKS  keyword. 

A  single  node  can  be  shared  by  any  number  of  links,  being  the  transmitter 
for  some  and  the  receiver  for  others. 

The  $NODES  keyword  line  can  be  followed  by  any  number  of  data  line 
pairs  which  specify  the  location  (latitude,  longitude,  altitude)  and  attri¬ 
butes  of  each  node.  Besides  its  location,  a  node  has  a  name,  a  number  and 
four  other  parameters:  transmitter  power,  receiver  noise  temperature,  re¬ 
ceiver  man-made  noise  class,  and  receiver  bandwidth.  If  a  node  is  used  as 
the  transmitter  for  a  given  link  then  the  receiver  parameters  are  ignored 
(and  vice-versa) . 

About  the  only  restriction  concerning  the  location  of  nodes  is 
that  they  must  not  be  too  near  (one  degree)  to  a  pole.  This  is  to  keep 


i 

* 

4 
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Table  6.  $NODES  keyword. 


Keyword  line: 
Data  Lines: 


$N0DES 

(A)  NODE,  LAT,  LON,  ALT,  TPOWER,  RNTEMP,  RMMNCL,  RBANDW 

(B)  HNODE 


Data  name 


Description 


NODE 

LAT 

LON 

ALT 

TPOWER 

RNTEMP 

RMMNCL 


RBANDW 

HNODE 


Node  number  [1,  1000] 

Node  latitude  (°N)  [-89,  89] 

Node  longitude  (°E)  [-360,  360] 

Node  altitude  (km)  [0,  20] 

Transmitter  power  (kw)  [0,  1000] 

0  =  use  TPOWER  =  1  kw 
Receiver  noise  temperature  (°K)  [0,  1000] 
0  =  use  RNTEMP  =  288  °K 
Receiver  man-made  noise  class  [0,  4] 

0  =  use  RMMNCL  =  2  (rural) 

1  =  remote 

2  =  rural 

3  =  residential 

4  =  industrial 

Receiver  bandwidth  (kHz)  [0,  1000] 

0  =  use  RBANDW  =  1  kHz 
Node  name  (24  characters,  6A4) 
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Table  6.  (cont.) 


Example 


SNODES  EXAMPLE 


SETUP  4  NODES:  SANTA  BARBARA.  BIG  FORK.  CAMBRIDGE  BAY.  AND  WASHINGTON. 


SNODES 


jNODE  lat 

LON 

ALT 

TPOWER 

RNTEMP 

RMMNCL 

RBAHDU 

1  34.25 

-119.41 

0 

l 

288 

2 

1 

SANTA  BARBARA. 

CALIF. 

2  48 

-114 

0 

1 

288 

2 

1 

BIGFORK,  MONTANA 

3  69 

-105 

0 

1 

288 

1 

1 

CAMBRIDGE  BAY, 

CANADA 

4  38.55 

-77 

0 

1 

288 

3 

1 

WASHINGTON,  D.C. 


/  HNODE 


Notes:  (1)  $NODE$  data  lines  A  and  B  are  paired 

(2)  Any  number  of  A/B  data  line  pairs  may  appear 

(3)  Node  numbers  must  be  unique 

(4)  Nodes  must  be  at  least  1  degree  from  any  pole 
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certain  geometric  approximations  from  blowing  up.  It  is  also  somewhat 
dangerous  to  define  a  link  which  crosses  directly  over  a  pole,  although  the 
code  does  not  check  for  this  occurrence. 

2.7  SLINKS  KEYWORD 

The  $LINKS  keyword  is  used  to  specify  the  attributes  of  any  number 
of  links.  Each  link  is  described  on  a  single  data  line  containing  six  para¬ 
meters:  link  number,  transmitter  node  number,  receiver  node  number,  iono¬ 
spheric  model  selector,  surface  type,  and  output  flag. 

Obviously  the  link  number  must  be  unique  and  the  transmitter  and 
receiver  nodes  must  be  defined  elsewhere  in  the  input  using  the  $NODES  key¬ 
word.  The  ionospheric  model  selected  is  up  to  the  users'  preference.  Pros 
and  cons  on  the  different  models  [Aerospace,  RADC  (polar),  and  ITS]  are 
given  in  a  companion  report  (Reference  2) .  The  surface  type  (land  or  sea) 
is  used  to  compute  ground  reflection  losses  on  multihop  paths.  It  should 
be  set  according  to  the  characteristics  at  the  midpoint  of  the  link.  The 
settings  provided  give  reasonable  results  in  the  absence  of  a  world  map  of 
surface  characteristics.  Finally,  the  print  flag  controls  whether  the  en¬ 
tire  ionospheric  data  table  (around  3  pages  worth)  is  printed  or  not. 

2.8  $HFCALC  KEYWORD 

The  major  calculation  performed  by  HFNET  is  the  HF  propagation 
calculation.  The  standard  HF  calculation  event  consists  of  several  related 
calculations  which  can  be  expressed  in  pseudo-code  as  follows: 

•  loop  over  links; 

•  initialize  for  this  link; 

•  loop  over  hops; 

•  compute  N-hop  MUF  and  decile  frequencies; 

•  compute  propagation  geometry,  losses,  signal  strength, 
noise  and  S/N  for  MUF; 
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Table  7.  $LINKS  keyword. 

Keyword  line:  SLINKS 

Data  line:  LINK,  TXNODE ,  RXNODE,  IMODEL,  SURTYP,  PFLAG 


Data  name 


Description 


LINK  Link  number  [1,  1000] 

TXNODE  Transmitter  node  number  [1,  1000] 

RXNODE  Receiver  node  number  [1,  1000] 

IMODEL  Ionospheric  model  [0,  3] 

0  =  use  IMODEL  =  1  (Aerospace) 

1  =  Aerospace 

2  =  RADC  (POLAR) 

3  =  ITS 

SURTYP  Surface  type  [0,  2] 

0  =  use  SURTYP  =  1  (land) 

1  =  land 

2  =  sea 

PFLAG  Print  flag  [0,  1] 

0  =  off 

1  =  on,  print  full  ionospheric 
data  table 
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Table  7.  (cont.) 


Example 


-SLINKS  EXAMPLE 


SETUP  3  LINKS:  (1)  SANTA  BARBARA  TO  BIG  FORK.  AEROSPACE 

(2)  UASHINGTON  TO  CAMBRIDGE  BAY.  RADC  (POLAR) 

(3)  SANTA  BARBARA  TO  UASHINGTON.  ITS 


SLINKS 


LINK  TXNQDE  RXNODE  IMODEL  SURTYP  PFLAG 


112  110 

4  3  2  1  0 

14  3  10 


Notes:  (1)  Any  number  of  SLINKS  data  lines  may  appear 


(2)  Link  numbers  must  be  unique 


(3)  Transmitter  and  receiver  node  numbers  must  correspond  to 
nodes  set  up  using  the  $N0DES  keyword 
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Table  8.  $HFCALC  keyword. 


Keyword  line: 
Data  line: 


$HFCALC 

Tl,  T2 ,  DT,  LINK,  FREQ,  HOPS,  EFLAG,  PFLAG 


Data  name 


Description 


Tl 

T2 

DT 

LINK 

FREQ 


HOPS 


EFLAG 


PFLAG 


First  event  time  (DHMS)  [0,  +] 

Last  event  time  (DHMS)  [0,  +] 

Event  rescheduling  interval  (DHMS)  [0,  +] 

Link  number  [0,  1000] 

0  =  use  all  links 

Frequency  (MHz)  [-30,  -2]  U  [0]  U  [2,  30] 
-  =  use  ABS  (FREQ)  only 
0  =  use  MUF  only 
+  =  use  MUF  and  FREQ 

Number  of  hops  (packed,  H^  H^)  [00,  66] 

00  =  use  internally  computed 
defaults  (see  notes) 

H-|  H2  =  use  H-|  to  H2  hops  (H^H2) 

Event  flag  [0,  1] 

0  =  normal 
1  =  multimode 

Print  flag  [0,  1] 

0  =  off 

1  =  on,  print  detailed  output 
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Table  8.  (cont. ) 


Example 


SHFCALC  EXAMPLE 


HF  PROPAGATION  CALCULATIONS  EVERY  HALF  HOUR  FROM  19:00  TO  21:00  i'GMT) 
AND  EVERY  HOUR  FROM  21:00  TO  19:00  (THE  NEXT  DAY). 

ALL  LINKS  --  FREQ  =  MUF.  13.  26  MHZ  --  DEFAULT  HOPS 
DETAILED  PRINTOUT  FOP  20  MHZ  ONLY. 


SHF CALC 


:  T 1 

T2 

DT 

LINK 

FREQ 

HOPS 

EFLAG 

PFLAG 

0  1  900 

02100 

0030 

0 

10 

00 

0 

0 

01900 

02100 

0030 

0 

-20 

00 

0 

1 

02200 

1  1900 

0100 

0 

10 

00 

0 

0 

02200 

1  1900 

0  1  00 

0 

-20 

00 

0 

1 

Notes:  (1)  Any  number  of  $HFCAtC  data  lines  may  appear 


(2)  Unless  zero,  the  link  number  must  correspond  to  a  link  set 


up  using  the  $IINKS  keyword 

(3)  The  defaults  for  hops  are:  1  hop  per  3000  km  (minimum)  and 
1  hop  per  1000  km  (maximum),  rounded  to  nearest  hop;  this 
results  in  the  following  table: 


Kilometers 

1000 

2000 

3000 

4000 

5000 

6000 

9000 

12000 

15000 

18000 


Min,  hops 

1 

1 

1 

1 

2 

2 

3 

4 

5 

6 


Max,  hops 

1 

2 

3 

4 

5 

6 
6 
6 
6 
6 
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•  compute  the  probability  of  mode  occurrence  for  speci¬ 
fied  frequency,  f; 

•  compute  propagation  geometry,  losses,  signal  strength, 
noise  and  S/N  for  f; 

•  next  hops ; 

•  next  link; 

The  $HFCALC  keyword  is  used  to  set  up  any  number  of  HF  calculation 
events.  Each  $HFCALC  data  line  contains  eight  parameters:  first  event  time, 
last  event  time,  event  rescheduling  interval,  link  number,  frequency,  mini¬ 
mum  and  maximum  number  of  hops,  event  flag  and  print  flag. 

The  "usual"  setting  for  the  link  number  is  zero.  This  causes  the 
calculation  to  be  performed  for  all  links  simultaneously.  This  feature  is 
very  convenient  when  there  are  many  links  in  the  problem.  However,  the 
flexibility  to  set  up  HF  propagation  calculations  with  different  times, 
frequencies  and  hops  for  each  link  is  available  if  required. 

The  frequency  parameter  is  interpreted  in  a  slightly  complicated 
fashion.  Since  the  MUF  is  a  very  important  parameter  in  most  HF  communica¬ 
tions  studies,  the  calculation  of  the  MUF  is  included  as  part  of  the  "normal" 
HF  calculation  event  unless  explicitly  suppressed.  The  MUF  calculation  is 
made  whenever  the  specified  frequency  is  positive  or  zero.  A  negative  fre¬ 
quency  suppresses  the  MUF  calculation  (the  absolute  value  of  the  frequency 
is  then  used),  while  a  zero  frequency  causes  the  MUF  alone  to  be  used.  If 
the  MUF  and  several  fixed  frequencies  are  desired  then  these  frequencies 
should  be  specified  on  separate  $HFCALC  data  lines  with  the  first  one 
positive  and  the  rest  negative  (to  suppress  the  MUF;  see  example  in  Table  8.) 

The  number  of  hops  parameter  is  somewhat  of  a  nuisance.  The 
recommended  value  is  zero,  which  lets  the  code  decide  how  many  hops  are 
feasible  depending  on  the  length  of  the  link.  However,  if  only  1  hop  paths 
are  of  interest  then  this  parameter  allows  the  user  the  option  of  suppressing 


multiple  hop  mode  calculations.  The  factors  used  to  internally  compute  the 
number  of  hops  to  try  are:  1  hop  per  3000  km  (minimum)  and  1  hop  per  1000  km 
(maximum),  rounded  to  the  nearest  hop.  These  factors  result  in  the  table 
shown  at  the  bottom  of  Table  8.  Note  that  the  hops  parameter  is  a  two 
digit  packed  number.  For  example,  HOPS  =  36  is  interpreted  as  "try  3  to  6 
hops". 


The  event  flag  is  used  to  select  between  the  "normal"  and  the 
"multimode"  HF  propagation  calculation.  Basically  the  multimode  calcula¬ 
tion  is  a  more  detailed  calculation  which  attempts  to  study  the  possibility 
of  multipath  propagation  by  1)  using  several  frequencies  (in  place  of  a 
single  MUF)  chosen  for  their  likelihood  of  producing  multipath  and  2) 
using  a  modified  propagation  algorithm  which  allows  more  than  one  ray  to 
propagate  at  a  given  frequency.  Except  for  these  differences  the  multi- 
mode  variation  is  very  similar  to  the  normal  HF  calculation.  Use  of  the 
multimode  option  is  recommended  only  for  single  links  and  for  no  more 
than  3  hops. 


The  print  flag  is  used  to  control  the  amount  of  detailed 
output  from  the  propagation  calculations.  When  the  print  flag  is  turned 
"on"  fairly  elaborate  output  is  generated  for  each  successful  mode  and  an 
explanation  is  offered  for  each  unsuccessful  mode.  Turning  the  print  flag 
"off"  still  leaves  the  summary  output  which  is  adequate  for  many  applica¬ 
tions.  See  Section  3  for  examples  of  both  the  detailed  and  summary  out¬ 
puts. 

2.9  $PLCALC  KEYWORD 

In  addition  to  conventional  ionospheric  skywave  propagation, 
HFNET  is  capable  of  computing  nuclear  induced  "plume  mode"  propagation. 

This  type  of  anomalous  propagation  involves  bouncing  HF/VHF  rays  off  of 
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Table  9.  $PLCALC  keyword. 


Keyword  line:  $PLCALC 

Data  line:  T1 ,  T2,  DT,  LINK,  FREQ,  PFLAG 


Data  name 

T1 

T2 

DT 

LINK 

FREQ 

PFLAG 


Description 

First  event  time  (DHMS)  [0,  +] 

Last  event  time  (DHMS)  [0,  +] 

Event  rescheduling  interval  (DHMS)  [0,  +] 

Link  number  [0,  1000] 

0  =  use  all  links 

Frequency  (MHz)  [0]  U  [20,  100] 

0  =  use  "MUF"  only 
+  =  use  FREQ  only 

Print  flag  [0,  1] 

0  =  off 

1  =  on,  print  detailed  output 


Example 


;  SPLCALC  EXAMPLE 
.■  '  —  —  —  —  —  — 

:  PLUME  MODE  CALCULATIONS  EVERY  HALF  HOUR  FROM  19:00  TO  21:00  (GMT) 
;  AND  EVERY  HOUR  FROM  21:00  TO  19:00  (THE  NEXT  DAY). 

:  ALL  LINKS  —  FREQ  =  "MUF"  —  DETAILED  PRINTOUT 


TPLCALC 


:  T 1 

T2 

DT 

LINK 

FREQ 

PFLAG 

01900 

02100 

0030 

0 

0 

1 

02200 

1 1900 

0  1  00 

0 

0 

1 

Notes: 

(1) 

Any  number  of  $PLCALC 

data  lines  may  appear 

(2) 

Unless 

zero,  the 

link 

number  must  correspond  to  a  link  set 

up  using  the  $LINKS  keyword 
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regions  of  enhanced  electron  density  created  by  nuclear  explosions.  Three 
types  of  plume  modes  are  computed: 

•  specular  reflection  off  of  the  side  of  a  plume 

•  isotropic  reflection  off  or  fireballs  and  plume  bases 

•  forward  scatter  off  of  plume  bases. 

Because  plume  mode  propagation  is  not  as  well  understood  as  normal  ionospheric 
propagation  and  because  it  is  highly  coupled  to  the  nuclear  phenomenology 
models,  the  results  of  plume  mode  calculations  need  to  be  analyzed  closely. 

The  models  should  be  considered  as  preliminary  versions  which,  while  useful 
when  used  wisely,  are  not  the  final  word  in  plume  mode  modeling. 

Plume  mode  calculations  are  organized  and  set  up  very  much  like 
HF  propagation  calculations.  The  $PLCALC  keyword  is  used  to  specify  any 
number  of  plume  mode  calculation  events.  Each  $PLCALC  data  line  contains 
the  following  parameters:  first  event  time,  last  event  time,  event  resched¬ 
uling  interval,  link  number,  frequency  and  print  flag.  Most  of  these  para¬ 
meters  act  just  like  the  corresponding  parameters  on  the  $HFCALC  data  line. 

The  frequency  parameter  is,  however,  interpreted  in  a  slightly  dif¬ 
ferent  way  than  the  $HFCALC  frequency.  Since  plume  modes  are  primarily  a  VHF 
phenomena  the  frequency  must  be  zero  or  between  20  and  100  megahertz;  nega¬ 
tive  frequencies  are  not  allowed.  If  entered  as  zero  the  code  will  compute 
the  "best"  frequency  to  use  depending  on  the  type  of  mode,  the  electron 
density  in  the  cloud,  the  available  scattering  cross  section  and  the  mode 
geometry.  If  the  frequency  is  positive  then  that  frequency  alone  will  be 
used. 

2.10  $DEBUG  KEYWORD 

HFNET,  like  any  large  computer  code,  may  contain  bugs.  Be  they 
incorrect  physical  models,  inappropriate  computational  algorithms  or  plain 
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Table  10.  $DEBUG  keyword. 


Keyword  line:  $DEBUG 
Data  line:  HDEBUG 

Description 

Names  of  up  to  10  subroutines 
(8  characters  each,  10A8) 

Example 

;  SDEBUG  EXAMPLE 

;  DEBUG  OUTPUT  FROM  MOD10  AND  PLMODE  (PLUME  MODE  SUBROUTINES). 

SDEBUG 
: HDEBUG 

MOD  10. . .PLMODE. . 

Notes:  (1)  At  most  one  $DEBUG  data  line  may  appear 

(2)  Subroutine  names  must  be  left  justified  in  fields  of  8 
characters  and  right  filled  with  periods 


Data  name 
HDEBUG 
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coding  errors  the  result  is  the  same:  the  code  may  not  work  correctly  for 
some  combinations  of  input  data.  Since  the  task  of  totally  eradicating  all 
bugs  in  any  large  code  is  clearly  impractical,  the  real  problem  is  to  make 
it  less  likely  for  a  bug  to  go  undetected  and,  once  detected,  to  make  it 
easy  to  locate  the  exact  nature  of  the  bug. 

The  detection  of  bugs,  especially  subtle  bugs,  is  a  difficult 
problem.  The  code  itself  trys  to  locate  possible  bugs  by  making  consistency 
checks  on  the  data  it  uses.  These  checks  are  of  the  grossest  nature  but  do 
occasionally  find  bugs  which  would  probably  have  gone  unnoticed  otherwise. 
However,  the  best  line  of  defense  against  bugs  (and  the  incorrect  results 
they  foster)  is  a  skeptical  user  of  the  code.  The  user  should  never  rely  on 
the  results  of  any  computer  code  without  first  checking  to  see  if  the  results 
are  reasonable,  and  if  they  are  not,  finding  out  why  not. 

Once  a  bug  has  made  its  presence  known  by  causing  the  program  to 
abort  or  by  giving  unreasonable  results,  the  problem  still  remains  to  locate 
and  fix  the  bug.  The  traditional  method  of  debugging  a  large  FORTRAN  code 
is  to  scatter  diagnostic  print  statements  throughout  those  parts  of  the  code 
where  the  bug  is  thought  to  reside.  Once  the  bug  has  been  eradicated  these 
print  statements  are  then  removed,  only  to  reappear  again  in  the  future  when 
some  new  bug  crops  up. 

HFNET  extends  this  traditional  debugging  method  by  1)  incorpora¬ 
ting  debug  print  statements  before  the  occurrence  of  a  bug  and  2)  leaving 
them  in  the  code  in  such  a  way  that  they  can  be  turned  on  or  off  easily 
(without  recompiling  or  relinking  the  program).  Because  of  this,  HFNET 
debug  output  can  be  used  not  only  to  find  bugs  but  also  to  provide  very 
detailed  output  which  the  regular  output  cannot  provide. 


The  $DEBUG  keyword  is  used  to  request  that  detailed  debug  output 
be  produced.  This  output  is  turned  "on"  or  "off"  on  a  subroutine  by  subrou¬ 
tine  basis,  with  up  to  10  subroutines  turned  "on"  at  the  same  time.  For 
the  most  part  this  output  is  only  useful  to  someone  reasonably  familiar  with 
the  code  and  normally  requires  a  source  listing  for  proper  interpretation. 
Table  10  shows  the  format  required  to  request  debug  output.  One  warning: 
debug  output  can  be  quite  lengthy,  so  it  should  be  used  with  caution. 

2.11  $RUN  KEYWORD 

The  $RUN  keyword  is  the  last  of  the  input  keywords.  It  is  the 
only  keyword  that  is  required  and  it  must  be  the  last  keyword  appearing 
in  the  input  file.  The  $RUN  keyword  has  two  functions:  1)  it  signals 
the  end  of  the  input  phase  and  the  beginning  of  the  simulation  phase  and 
2)  it  allows  the  user  to  specify  when  the  simulation  should  stop.  After 
the  $RUN  keyword  has  been  processed  the  input  module  returns  to  the  simula¬ 
tion  manager.  If  no  errors  were  detected  in  the  input,  the  simulation 
manager  will  then  proceed  to  execute  events  until  the  stop  simulation  event 
is  encountered. 

The  $RUN  data  line  contains  two  parameters:  the  simulation  stop 
time  and  an  output  flag.  The  output  flag  is  used  to  control  printouts  of 
the  event  list  and  dynamic  storage  at  the  beginning  and  end  of  the  simula¬ 
tion. 
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Table  11.  $RUN  keyword. 


Keyword  line:  $RUN 
Data  line:  TSTOP,  PFLAG 


Data  name 

TSTOP 

PFLAG 


Description 

Simulation  stop  time  (DHMS)  [0,  +] 

Print  flag  [0,  1] 

0  =  off 

1  =  on,  print  event  list  and 
dynamic  storage  maps 


Example 


;  JRUN  EXAMPLE 

;  STOP  SIMULATION  AT  19:01  (GMT)  THE  NEXT  DAY;  PRINT  EVENT  LIST  AND  D.S.  MAPS. 
$RUN 

; TSTOP  PFLAG 

11901  1 


Notes : 


(1)  The  $RUN  keyword  must  appear  and  be  the  last  keyword  in 
the  input  file 

(2)  Exactly  one  $RUN  data  line  must  appear 
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SECTION  3 
HFNET  OUTPUT 


3.1  GENERAL 


Printed  output  from  HFNET  can  be  divided  into  two 
1)  detailed  simulation-time-ordered  output,  and  2)  summary 
output.  This  section  describes  these  two  types  of  printed 
detail . 


categories : 
post-processed 
output  in 


There  are  three  basic  differences  between  the  detailed  output  and 
the  summary  output.  The  most  obvious  difference  is  that  the  summary  output 
contains  much  less  information  than  the  detailed  output.  For  example,  de¬ 
tailed  output  describing  a  single  HF  skywave  mode  includes  printouts  of 
raypath  geometry,  propagation  losses,  noise  sources  and  signal/noise.  All 
this  can  easily  amount  to  an  entire  page  of  output  depending  on  the  number 
of  hops.  In  contrast,  the  same  mode  is  described  on  a  single  line  in  the 
summary  output. 

The  second  difference  is  the  "time"  at  which  the  outputs  are 
created.  Simulation-time  output  is  generated  "on  the  fly"  during  the 
course  of  the  simulation.  Due  to  the  fact  that  events  are  executed  in  a 
strictly  time  ordered  sequence,  printouts  from  many  different  events  can 
be  interleaved.  Consequently,  a  good  deal  of  "page  flipping"  through  the 
output  to  find  the  printout  of  interest  is  commonplace.  On  the  other  hand, 
post -processed  output  is  generated  by  a  separate  program,  HFPOST,  after  the 
main  program,  HFNET,  has  finished  execution.  This  enables  the  summary  output 
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to  be  sorted  and  reformatted  in  a  way  that  is  convenient  for  the  user  to 
read  instead  of  just  convenient  for  the  program  to  generate. 

The  third  difference  is  that,  for  the  most  part,  the  detailed 
simulation-time  output  can  be  suppressed  using  the  various  print  flags 
described  in  Section  2.  The  summary  output,  on  the  other  hand,  is  automa¬ 
tically  generated.  (Of  course  the  user  can  elect  not  to  run  the  post¬ 
processor  program  and  thereby  suppress  all  summary  output.)  The  output 
from  a  "typical"  HFNET  run  consists  of  some  user  selected  detailed  output 
in  addition  to  the  summary  output. 

As  a  means  of  conveniently  discussing  HFNET  printed  output,  exam¬ 
ples  will  be  taken  from  a  test  problem  that  shows  most  of  the  major  types 
of  output  as  well  as  exercises  a  large  part  of  the  code.  A  listing  of  the 
input  file  for  the  test  problem  is  shown  in  Table  12.  Except  for  $DEBUG, 
all  of  the  keywords  discussed  in  Section  2  are  used  in  this  input  file. 

The  test  problem  is  a  calculation  of  nuclear  effects  on  three 
specific  links  shortly  after  the  detonation  of  two  nuclear  devices.  The 
two  nuclear  bursts  are  a  high  altitude  megaton  class  burst  over  the  western 
United  States  followed  closely  in  time  by  a  low  altitude,  large  yield  burst 
over  the  eastern  part  of  the  country.  The  links  in  the  problem  are  1)  Santa 
Barbara,  California  to  Big  Fork,  Montana,  2)  Washington,  D.C.  to  Cambridge 
Bay,  Canada  and  3)  Santa  Barbara  to  Washington.  The  test  problem  includes 
both  HF  propagation  and  plume  mode  calculations  made  approximately  30  and 
90  minutes  after  the  burst.  Figure  1  shows  the  problem  geometry  including 
the  location  of  each  of  the  two  bursts,  three  links  and  four  nodes.  A  pro¬ 
spective  user  should  at  this  point  study  Table  1  carefully  (reviewing 
Section  2  if  needed)  until  he  or  she  fully  understands  how  each  keyword, 
data  and  comment  line  is  being  used. 


36 


Table  12.  Test  problem  input  file 


:•  USERMAN.DAT  • 


SI  PENT 

; PROBLEM  IDENTIC [CATION 


example  problem  for  hfnet  USER'S  manual 


same  I ENT 


:  YEA* 

2  9?  4 

MONTH 

4 

DA'' 

] 

SSN 

113 

*e 

4 

SOLAR 

VED1ND 

0 

S/NUJND 

0 

WINDS 

1 

seursts 

:TIME 

LA* 

uON 

ALT 

Y 

PY 

XY 

spread 

CBFLAG 

PBFLAG 

o.esi 

40 

-120 

200 

1 

C 

0 

c 

0 

16 

t:zt2 

uo 

-80 

1 

20 

0 

c 

0 

0 

0P 

SCLDOJT 


T*  T< 

ei?ee  11900 

DT 

0I0C 

PBFLAG 

1C 

SNQPES 

;  NOPE  LAT 

LON 

ALT  TPOWER 

RNTErf» 

PWCl  ReANDU  /  NODE  NfirE 

34.25 

0  1 

282 

3  J 

SmnTA  EAR&APh, 
42 

.  CAL  If- 

-114 

0  1 

ZE8 

2  1 

2 ItfPORK .  MONTANA 

3  €9  -105 

0 

288 

1  ! 

CAM&FIDGE  BAY- 
4 

.  CANADA 

0  1 

2PS 

4  1 

OlASHlMGTCl!-.  t. 

,  c . 

XL  INK; 


s^|N»  TYNCDE 

RXNODE 

I MODEL 

5‘JRTVP 

pPL  AG 

1  1 

2 

j 

1 

C 

2  4 

T  ' 

L 

J 

1 

3  1 

*• 

' 

0 

SHF CALC 

:T1  T2 

vr 

l.  INK 

FRSC 

HOPS 

EeLAG  pflag 

Cl 906  11900 

0100 

C 

10 

0C 

C  0 

81900  1 1900 

6100 

6 

-20 

60 

e  1 

splCAuC 

;T1  T2 

I'T 

LINK 

FPEO 

pflag 

01900  11906 

0100 

0 

0 

1 

SPUN 

:TSTOF 

PPI.H& 

C2C01 

1 
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3.2 


SIMULATION  START  OUTPUT 


The  initial  output  from  the  HFNET  code  is  shown  in  Table  13. 

This  output  begins  with  a  short  header  containing  the  version  of  the  code 
being  used,  the  name  of  the  input  file  and  the  date  and  time  of  the  run. 
Following  this  is  a  listing  of  the  entire  input  file  enclosed  in  a  box  of 
asterisks.  If  any  input  errors  were  detected  then  error  messages  would 
also  appear  interspersed  with  the  echoed  input.  Finally,  printouts  showing 
the  $AMBIENT  and  $DEBUG  data  being  used  during  this  run  follow  the  echoed 
input. 


Table  14  shows  the  event  list  and  dynamic  storage  printouts  at 
the  beginning  of  the  run  (after  all  input  has  been  read  and  stored  but  before 
any  further  events  have  been  executed).  These  printouts  are  optional  and 
appear  here  because  the  print  flag  on  the  $RUN  data  line  was  turned  "on". 

The  headings  for  the  event  list  printout  are: 


EVENT 

-  event 

sequence  number 

TIME 

-  event 

simulation  time  (GMT  in  DHMS  format) 

TYPE 

-  event 

type  number  (see  Appendix  A) 

NAME 

-  event 

descriptive  name 

EVENT  DATA  -  16  words  of  event  data  for  each  event 

Notice  that  the  first  four  events  all  have  a  simulation  time  of 
-0.01  (-1  second).  These  events  are  initialization  events,  they  read  in  the 
wind  model  database  and  compute  ionospheric  data  tables  for  the  three  links 
in  the  problem.  Following  the  initialization  events  are  seven  events  which 
were  explicitly  set  up  using  the  $BURSTS,  $CLD0UT,  $HFCALC,  $PLCALC  and  $RUN 
keywords  in  the  input  file. 
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Table  13.  Simulation  start  output. 


.*?* * ♦  *  ***.*  *  *  *  **#**'K*  ****;***  *  ***  tm*  *  *•  4-  lot  ***  ** 

*  HFtlET  *  VERSION  5  (VAX)  *'  SEPTEtflFP  1S79  * 

M’*  **  4-  f-T*  fir  |»f »  l,  *****  *  *#  M*  *  **  K**4  *«««*  4.4  f** 

SIMULATION  INPUT  FILE:  DATAsUSERMhN. DPT 
I*hTE  AMI*  TIME  OF  PUN:  31-0CT-79  17:10:16 


♦  :*  USEPMwN.DAT  » 

■M 

4  ; - - - -- - - - 

♦  5 1  LENT 

♦  ; - 

'  :PPOBlEM  IDENTIFICATION 


*  SAMPLE  PROBLEM  FOP  Hr  NET  USER'S  MANUAL 

4 

■♦  SAME  I ENT 

f  ;  T 'EhP  MONTH  LAY  SSN  KP  SOLAR  VEUINt*  VNUINL  MINI'S 

*  1524  4  1  113  4  D  0  0  3 


>. 

A 

S 

4- 

A 

*. 

A 

A 

A 


■4  ; - - - -  —  - — - - - — - - - - - - A 

♦  SBUPSTS  * 


4  : TIME  LhT  LON  ALT  *Y  MY  SPREAD  CBFLAG  PBFLAG  * 

*  81S31  40  -120  200  1  0  0  0  0  10  * 

*  01232  40  -SO  1  20  0  0  C  0  G0  ■* 


SELLOUT 
*  •  — - 

*  :Ti  T2  I*T  PBFLAG 

4  01900  11900  0100  10 


A- 

A 

A- 

A 

.4 

A 
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Table  13.  (cont.) 


■ 

SHODEi 

4- 

‘ 

;NODE  LhT 

-ON 

alt 

TPOWEP 

PNTFMP 

Prime  L 

RBANDU  /  NODE  NAME  * 

•4. 

1  7-4. 2S 

-i 19.41 

0 

1 

208 

3 

1  * 

* 

SANTA  BAFBAPA. 

Chi.  IF. 

» 

4 

2  48 

-114 

0 

1 

269 

* 

1  4 

* 

B1GF0PI  .  MONTANA 

4 

*• 

3  bS 

-105 

6 

1 

28S 

: 

4. 

4. 

CAMBRIDGE  BAY. 

CANADA 

*• 

4 

4  72 . 55 

-77 

ri 

1 

28? 

4 

1  * 

4 

LUSHING TOM.  0. 

c. 

4 

* 

4 

4 

iuINKS 

4 

4. 

> 

4 

: L I NK  T^NODE 

R-.-NOt’E 

I MODEL 

S11PTYP 

PFLAG 

i. 

4' 

1  1 

L 

1 

1 

o 

•4- 

4 

2  4 

2 

2 

1 

: 

* 

4. 

3  I 

4 

3 

i 

0 

4 

4. 

4' 

:♦ 

SHFCALC. 

> 

-f 

:T|  72 

PT 

LINK 

FREO 

HOPS 

EFL9G 

PPLAG  » 

* 

0 1 900  11900 

0100 

0 

10 

00 

A 

o  •* 

01900  11900 

0100 

0 

-20 

00 

0 

1  * 

4 

JPlChlC 

4 

:T!  ~2 

DT 

LINK 

FPEC 

PFLAG 

4 

* 

4* 

4 

01900  11909 

0100 

c 

0 

: 

TPM*. 

4 

; TCTO- 

pflag 

4 

0200 1 

1 

/ 

*-4c*o*r***4c«c  ryy* 4o*r**c**>ry*:H:*^^**4*f  a 

**  AMBIENT  LEATHER  CONDITIONS  ** 


’v'EAP 

1084. 


MONTH 

4. 


DAY 

1. 


SSN 

IH. 


KP 

4. 


SOLAR 

0. 


VEWIND 

0.06 


VNWIND 

0.00 


**  DEBUG  OUTPUT  HAS  BEEN  REQUESTED  FROM  THE  FOLLOWING  SUBROUTINES  ** 
NONE .... 


WINDS 

t. 
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Table  14.  Event  list  and  dynamic  storage  printouts 


*  HFNET  EVENT  L  1ST  CnHTulNS  tl  CVTIITS  *  THE  f  1PC-  f  II  F  VTHTv.  APR  LISTED  BEL OM  * 


EVENT 

Tilt 

TYRF 

NAME 

EVENT  Dof  A 

i  t) 

-iT.ftl 

IS 

RFhDUIND 

-  1  .  OOOF  All  I 

l.SOOEnll 

0.0001  tfllj 

o.onot  +on 

O.AflOE+00 

0.0G0E+00 

0.B00E+00 

0.000E+O0 

H.H+TUl  <00 

tl .  +JOUt  ♦  on 

H.  01101  +00 

o.ofiuE  +00 

0. AA0F+0O 

0  .  LluAE+UU 

O . OOOE +00 

0 . OoOE +mm 

•  2) 

-0.01 

iHNii  PHR 

1  .  iJi u ip  <00 

'.‘.WOOL  +00 

l  .Oi  nit  +00 

0 .W00E +00 

+.1 .0001  +00 

0  .Of IMF.  +  00 

0. 0006+00 

0 . 00OE  +00 

0.  Mill ll.  +00 

0 .  CniOE  *  wfl 

0.0H0E tOW 

0 . 0+JUL  +00 

0.00 HE +00 

0 .  tU+WE+UU 

0 . 0OUE+0U 

U - OOOE  +00 

1) 

-  0 .  A  l 

? 

liJHO.F-Hk 

-  1 .UUUF »00 

r.UUOi  i AO 

.‘.(tiiOE  ttiU 

1 . 0uer i 00 

0.00!  IL0 10 

0  .tmuC+oo 

0  .  ui  10C  +00 

0 . U0OL +  00 

o. m-titL' » ufj 

ti.  Atif.E  +  00 

U.OHWl  +00 

0.000F  *0u 

U  .  tJObE  + 110 

O.0Li0E+fJU 

0.00+TE  +  OO 

0 . WOOL +0W 

i  <r> 

-0.01 

1  UNO  4' HR 

* l .HOUR  too 

r.umiL +uo 

3.O0UE  * (10 

o.oouf  mu 

O.00OE  +00 

O.00OE+0O 

0 . 00 WE +00 

0.0+101  +00 

O.HHHE  *Ut) 

O.000E  »llO 

0 . OO0E  +H0 

u.niruc  +00 

0 .  CiHiTE  +  00 
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In  the  dynamic  storage  memory  map  only  three  data  lists  are  shown 
the  EVNT  list  containing  11  events?  the  NODE  list  with  datasets  for  4  nodes 
and  the  LINK  list  containing  3  link  datasets.  The  "TOTAL  CALLS"  and  "CPU 
TIME"  printouts  show  how  many  calls  have  been  made  up  to  this  point  to 
dynamic  storage  routines  and  how  much  CPU  time  (in  seconds)  it  took  to 
execute  those  calls.  These  numbers  are  very  useful  in  evaluating  the 
amount  of  overhead  the  dynamic  storage  system  entails.  If  this  overhead 
is  too  great  then  some  alteration  of  the  dynamic  storage  parameters  is 
probably  needed. 

3.3  IONOSPHERIC  MODEL  OUTPUT 

The  complete  ionospheric  data  table  computed  for  the  Washington 
to  Cambridge  Bay  link  using  the  RADC  (polar)  model  is  shown  in  Table  15. 

The  headings  in  the  top  portion  of  this  table  are  self  explanatory.  The 
body  of  the  table  shows  ionospheric  parameters  above  seven  equally  spaced 
points  on  the  great  circle  path  from  the  transmitter  to  the  receiver.  At 
each  of  these  points  data  are  given  for  three  ionospheric  layers  (E,  FI  and 
F2)  and  for  25  times  during  the  day.  The  headings  for  this  table  are: 

PLACE  -  place  number  (1-7);  latitude  (°N)  and  longitude  (°E) 
of  this  place  (place  1  is  the  transmitter;  place  7 
the  receiver) 

TIME  -  local  time  (hours) 

ELECTRONS  -  layer  peak  electron  density  (cm  ) 

C-FREQ  -  layer  critical  frequency  (MHz) 

ALT  -  layer  altitude  (km) 

1/2  THICK  -  layer  semi-thickness  (km) 

Ionospheric  data  tables  like  the  one  shown  in  Table  15  are  computed  and 
stored  for  every  link  set  up  using  the  $LINKS  keyword.  The  printing  of 
the  table  is  optional  and  controlled  by  the  $LINKS  print  flag. 

*• 

To  save  space,  only  the  first  10  datasets  in  any  dynamic  storage  list 
are  shown. 


Table  15.  Ionospheric  data  table  output 


*  •  f«r«ci  .*■*>>  ttyt:*  *■'.?&  «-(c  :’i 

*'  IONOSPHERIC  IhTh  r^*Bi_E  CPEmT£I>  ecfc 

*  *  ..«r.!  »  V  t.K-  “f  *  f 

1  *•*;('  **  ►..yy*' '*•*•*'>'*.*'  *•***"*  v-nnuymjir  y  r  *»*  *'<«. 

•id  :-E  uCCATIQN 

LAT  L3M 

1  ?SG)  «DEG' 

HLT 

urn 

AZIMUTH 

•  DEG » 

RANGE 

*  KM  1 

MODEL 

4.  WASHINGTON.  D.C. 

T3.55 

-77.00 

0.00 

3*42 . 4? 

3773.  ]5 

RQDC  ( POLAR V 

lAMBPIWE  Bh  •’.  CANADA 

59.00 

-105.0*3 

0.00 

1 “S.92 

3 : 

JEP 

Tt“F  gANDUlDTH 

MAN-MADE 

SL’PFACE 

'l;  * 

DEG  *> 

•HZ'. 

NO ! 

SE 

T*iPE 

2’?S .  00 

!  COO. 90 

REMO 

T6 

LmNI» 

6  L^'^EP 

y.15 m 

1  '2 

fc-r+y*- 

F !  LAYER 

TT'Kir* 

F2  LftMEP  tw»+ 

\s2 

31.  ACE 

TIME 

ElEC’PCMS 

C-r»£Q 

mLT 

THICK 

ELECTRONS 

C-FPEO 

alt 

THICK 

ELECTPON3 

C-FPEO 

ALT 

THICK 

LhT. 

•  HP 

CM-3' 

•.I1HZ, 

U.M> 

»  KHj 

•*Crf-3i 

'•MH2; 

U  f*? 

N1) 

' CM- 3 ' 

'MHZ> 

u:m> 

(r.n- 

0.00 

:.:?oe+04 

1.50 

115.31 

20 .  CO 

1 .0O0E+02 

9.03 

100.00 

40.00 

4.0155+05 

5.69 

367.35 

98.13 

;c.;: 

i .  ?e 

2 .  r?oe*t»4 

1.50 

i 15 . 31 

20.00 

1.0006+02 

0.0? 

100.00 

40 .  OO 

3.622E+05 

5.48 

376.27 

95.82 

- "7 . oo 

2 . 00 

2 . "30c +04 

1.50 

115.31 

20.09 

I.OO06+02 

O.03 

100.00 

40 . 00 

3.3796+05 

5.22 

376.37 

37 .  .32 

7.00 

: .  r  30E-H04 

! .  50 

115.31 

20.00 

1 . 0906+62 

0.0? 

100.00 

40.00 

3. 0316*05 

4.94 

378 . 37 

07.30 

J .  00 

2.790c +04 

1.50 

113.31 

20 . 00 

l . 0005+92 

O.09 

100.00 

40 . 00 

2.5165+05 

4.50 

>J<  6  .  Jr 

37.22 

5.00 

2 . “90E+04 

1.50 

115.31 

2C..30 

1 .0006+02 

0.03 

100.00 

40 . 00 

2. 3816+05 

4 . 33 

365.36 

36.25 

= .  00 

7.  -?15+04 

1 .  “5 

126.39 

20 . 00 

I .9006+92 

0 . 09 

10*3.00 

40 .  C0 

3 . 2446+05 

5.  11 

259.72 

64.51 

r .  oo 

l .  739E+05 

2.31 

1 14.63 

20.00 

1 .0006+02 

0.03 

100.00 

40 . 00 

5.241E+0S 

8.50 

236.22 

53.43 

■B .  00 

» .  36+05 

7.7? 

113. “7 

20.00 

1 . 00*36+62 

0.0? 

100.00 

40 . 00 

7.535E+05 

7 .  r? 

243.37 

61.13 

? .  00 

l . 65  4E-05 

3 . 65 

113.73 

20.00 

2.5476+95 

4. S3 

197.36 

40.00 

9. 1946+05 

3.61 

251.71 

63.32 

:o.oo 

1 .6326+05 

7.63 

114.0b 

20.09 

2.765E+05 

4. 72 

ISS.63 

40.00 

I . 045 E +06 

9.  13 

253.45 

66.6  4 

i  i .  00 

1 .77  IE +05 

1  -o 

! 14. *2 

20.O0 

2.513E+95 

4.35 

183.36 

40 .  CO 

1 . i “66+06 

9.74 

267.19 

69.34 

i : .  oo 

1.3205+05 

7.33 

1 14.52 

20.00 

2.9?:3c-05 

4.32 

131.31 

40.06 

1 . 282E+06 

10,  17 

274.94 

“1.95 

1 7 . 0*3 

1  •  5.006 +05 

3 . 33 

: 14.53 

2O.00 

3. 0276+05 

4.74 

132. 13 

40.00 

l . 3135+06 

10.2? 

232.68 

74.3  3 

14.00 

l."3E+0'5 

3  .  “  : 

1  14.30 

20.00 

2.998E+05 

4 .  ?  1 

136.1  l 

40.00 

1.2SSE+06 

10.  I? 

290 . 42 

76.61 

15.  5C 

1 .“02E+05 

2.70 

I 14.20 

20.00 

2. 2*7?E+05 

4.82 

193.65 

40.00 

1.2576+06 

10.07 

298.  17 

73 .  '-7 

•  ^ .  .;<o 

.  656c +05 

3.65 

114.03 

29.00 

2. 637 E +95 

4 . 63 

206 . 3 l 

40.00 

1 . 2295+06 

3 . 36 

305.31 

60.  14 

1 7 . 00 

1 . 6 15c +05 

3.61 

I  i  4 .  2  c 

20.O0 

2 . 7635+04 

1.49 

223.20 

40.00 

l. I79E+06 

3.75 

313.65 

91.33 

13.00 

1 . 3206+05 

3.26 

116.  15 

20.00 

1 . 0006+92 

0.O9 

100.00 

40.00 

l . 0716+05 

3 1 2? 

32 1 . 39 

100.56 

13.00 

2 . 7905+04 

;  .50 

1 15.21 

20.00 

1 .030E+92 

? .  09 

IOO.*?0 

40 . 00 

8.9556+05 

3' 50 

323.  14 

1 14. 14 

20 . 0O 

2.. “906+04 

1.50 

1  15.3! 

20.90 

1.0006+62 

0 . 09 

100.00 

40.00 

7 . 0066+05 

7.51 

326.33 

109.75 

21.00 

2 . 790E*04 

1.50 

1 15.71 

20.00 

t .OOOE+02 

0.09 

100.00 

40.00 

5.5646+05 

6.71 

344. 62 

104.40 

'  "* .  >0 

I. 7905+04 

1 . 30 

115.31 

20.09 

1 .0006-92 

0.O9 

100.00 

40.08 

4 . 5  7 1 E -05 

6. 2“ 

332 . 37 

101 .55 

2  2 . 00 

2.  “90E+04 

1.50 

115.31 

20 . 09 

1 . GO0E+O2 

9.09 

100.00 

40.00 

4. 463E+05 

6 .  OO 

360 .  1  1 

100.63 

24.00 

2 . 7905+04 

1.50 

115.71 

20.OO 

1 .0005+92 

9.09 

100.00 

40 . 00 

4. 0156+05 

5.63 

“67 . 65 

93.  13 

“i 

0 . 00 

2.  “C-05+04 

1.30 

1 15.21 

20 . 00 

1 .0906+02 

0.09 

100.00 

40 . 00 

3.315E+05 

5.  17 

369.51 

37.44 

42^3“ 

l.?0 

2. “505*04 

1  -7 

115.31 

29 . 08 

I . 0O0E+O2 

0.0? 

100.00 

40.00 

2.3796+05 

4.32 

453.10 

99.76 

-79.76 

:e 

2 .  “90E+04 

1.50 

115.31 

20.00 

1 .0006+02 

0.09 

100.00 

40 . 00 

2. 5596+05 

4.54 

453.  10 

99.61 

3 . 00 

2 .  “'306+04 

1 .  SO 

115.71 

20 . 00 

1 .OOOE+02 

9.09 

100.00 

40 . 00 

2.  1666+05 

4.20 

453. 10 

99.59 

4.00 

2. 7506+04 

1.50 

1 15.31 

20.00 

1 .OOOE+02 

0.69 

100.00 

40.O0 

l  .7S2E+Q5 

2.79 

453. 10 

99 . 60 

5.00 

2 . 73OE+04 

l.SO 

l 15.31 

20 . 00 

I .OOOE+02 

9.0? 

100.00 

40 . 00 

1 . 722E+05 

3.73 

453. 13 

99.71 

-: .  oo 

2.7906+04 

1.70 

115.3! 

20 . 00 

1 .O0O6+02 

0.O9 

100.00 

40.00 

2. 522E+05 

4.51 

275.67 

68.49 

"  .  00 

1 .4725+05 

3.44 

115.35 

20.00 

l.OOOE+02 

0.09 

100.00 

40.00 

3. 9536+05 

5.65 

239.46 

59.40 

3 . 00 

1.0175+05 

3.61 

l 13.93 

20 .  CO 

l . 0906  +02 

0 . 09 

100.00 

40 .  OO 

5 . 7 485+05 

6.8! 

“48.62 

61.59 

3 . 00 

1 .3305+05 

3.57 

114. 17 

20.00 

2.4156+05 

4.41 

202.74 

40.00 

7 . 254E+05 

7.65 

252.36 

64.  33 

10.00 

1 .626E+05 

3.62 

l 14.7? 

20 .  OO 

2.6416+05 

4.61 

192.7? 

40 . 00 

S. 4636+05 

8.26 

261.11 

86.?? 

.  1 .  *?c 

l .6?lE+?5 

2.63 

i 14.56 

20.00 

2.“?  66  +*35 

4.75 

137.05 

40 . 00 

9.613E+05 

9.30 

263.05 

63.65 

■  2 .  i'o 

’..“755-35 

3. “4 

1 14.65 

2O.0O 

2.399E+OS 

4.?3 

;?4. 55 

40 . 00 

1  .Q5“E+06 

?.23 

278.3? 

-7 

1  ~ .  00 

1.7435+05 

2.75 

1 14.65 

20.00 

2.52:6+05 

4. 65 

134.96 

40.00 

1  .  103E+O6 

3.45 

264.24 

74.86 

14.00 

1 .“136-05 

■?  -  *• 

1 14.55 

29.09 

2.595E+05 

4.33 

183. 33 

4*3 . 00 

1.  1216+08 

3.51 

292.08 

“6 . 38 

5  .  OO 

.  6665 +03 

3]  66 

114.30 

20.09 

2 .70 16+05 

4.75 

135 .  13 

40 . 00 

1 .  1156+06 

9.48 

299. S2 

73.  35 

1 6 .  CO 

1 .6565+05 

2.65 

3  1 4. 00 

20 . 08 

2 . 3J7E-05 

4.50 

205.59 

40.00 

1.096E+06 

9.40 

307 .57 

30.4? 

17.00 

: .  *045+05 

3.  “1 

1 14.06 

29.00 

5 . 23 15+04 

2.06 

225 . 95 

40.00 

1.0566+06 

9.23 

215.31 

91.7“ 

10.00 

1.5035+05 

3.4S 

115.7" 

20.00 

1.0006+92 

0.39 

100.00 

40.  CO 

9.  “306+05 

3.36 

323.05 

92.12 

;  O  .  ,’;j 

2.  "  306-04 

1.50 

1 15.21 

20.09 

l  .0*306+02 

3.09 

100.00 

40.00 

3.3236+05 

3.1? 

333. “9 

112.6? 

20 . 00 

2 .  “'-0E+O4 

1.39 

:  15.31 

20.00 

1 .0006+02 

0.0? 

100.00 

40.00 

6 . 6446 +05 

7. 32 

233.54 

105. 30 

21.00 

2 . 7 30t +04 

:  .50 

115.21 

29 . 00 

1 . 0*306+02 

0.O5 

100.00 

40 . 03 

5.278E+05 

6.52 

246.23 

i 02 . 66 

22.00 

2. "306+04 

:  .50 

115.31 

20.00 

1. OOOE+02 

0.0? 

100.00 

40.00 

4. 426E+05 

5. 37 

354 . 02 

37.  72 

2  2 .  OO 

2 . 7  = 06+04 

4  .SO 

1 15.21 

20.00 

1 . 9006+02 

0 . 03 

100.00 

40.00 

3.347E+05 

5.57 

361.77 

92.56 

24.00 

2 . "395+04 

1.50 

115.31 

20.90 

t .OOOE+02 

0.09 

100.00 

40 . 00 

3.315E+05 

c  •  f 

369.51 

97.44 

44 


rum* 


Table  15.  (cont.) 


< 


- 

0 . 00 

2.  :35E+*'5 

4. 15 

1 15.4? 

20.00 

1 .000E+O2 

0.09 

1O0.50 

40 . 00 

4.5566+05 

6 . 06 

•"0.36 

83 . 02 

4?.  25 

; .  oo 

. 2 01 £+05 

4.2: 

113.21 

20.00 

1 .0OOE+02 

5.09 

100.00 

48.00 

4.  1936+05 

5.82 

4Z3.3-: 

?7 . 4 1 

-32 . 1? 

: . 

2.  1 736*05 

4. 19 

115.2" 

20 . 00 

1 .00OE+02 

0.09 

100.00 

40.50 

3.3726*05 

5.5? 

433. ?3 

97.24 

7 .  oo 

'.  “"0£*04 

1.50 

115.31 

20.00 

1  .U0OE+02 

0.09 

100.00 

JO  .  00 

7. 3476*05 

5.1? 

433 . 93 

97.20 

4\  00 

.  730E+O4 

1.5C 

1 15.31 

20.00 

t . 090E+02 

0.09 

100.00 

45 . 00 

2 .5566 +*35 

4.54 

433.33 

97.21 

T.CC 

:.~?OE*04 

1.50 

1 15.71 

20 . 00 

l .  000E+92 

0 . 59 

100.00 

40.00 

1 .918E+05 

1  3T 

433 . ?8 

?"  .2? 

*5.00 

. 7*06+04 

1.50 

1 15.7 1 

20.00 

1 . 000E+O2 

0 . 09 

100.00 

40 .  '30 

2.0126*05 

4.03 

2?0 . 20 

- -*  37 

7 . 00 

.5466*05 

3.53 

l 15.60 

20.00 

1.O0OE-02 

0.09 

100.00 

40.00 

3.002E+G5 

4.32 

242.30 

60.41 

3 . 00 

.  593E-05 

3.53 

1 14.22 

20 . 00 

1 . 0O0E*02 

0.0? 

100.00 

40.00 

4. 3346*05 

5.?! 

246.37 

6  1 .  ?" 

3 .  OO 

.5656*05 

3.55 

t  14.12 

20.00 

1 .9006*02 

0.39 

100.00 

40.00 

5. 61  "6*05 

6. "3 

254.72 

04.65 

0.00 

.57bE*0S 

J.5e 

1 14.32 

20 . 00 

2. 499E+05 

4.49 

1*7.39 

40 . 00 

6.7226*05 

7.36 

262.46 

67. 26 

1 . 00 

.6196*05 

3.61 

*  14.49 

20.00 

2.6626+05 

4.63 

!?1 .55 

40 . 00 

7 . 7076*05 

7.38 

2 "8.20 

09.37 

2 .  CO 

. <5K£-05 

s  .85 

i 14.54 

20.00 

2. 7646+05 

4-72 

188.53 

40.00 

3.5516*05 

3.30 

2"  7 .  *3C 

72.41 

-  oo 

. Ob? E-OS 

3.07 

1 14.52 

20 .  CO 

2.5076+05 

4.76 

138 . 39 

40.00 

9.201E+05 

3.61 

285 . 6? 

74.32 

4.00 

.t'5rc+C5 

3. -35 

1  14.42 

20.00 

2.7916*05 

4.74 

191 . 17 

40 . 00 

9. 620E+C5 

3.21 

2? 3. 47 

"7.04 

5 . 00 

.0406*05 

3.04 

114.10 

20.00 

2 . 7126+05 

4.68 

197.  12 

40.00 

9.7506*05 

3.37 

301 . 13 

79.03 

a.  00 

.667E+05 

3.67 

113.53 

20.00 

2. 553E+05 

4.54 

207.32 

40.00 

?.603E*05 

3.30 

388.92 

80 . 75 

r .  00 

. 733E-05 

3 . 75 

113.79 

2O.0O 

7.4156*04 

2.44 

224.23 

45 . 00 

9.2676*05 

8.64 

3l0.6b 

32.08 

3.00 

.655E*C5 

3.65 

115. .39 

20.00 

* . O0CE+O2 

0 . 39 

IC0.O0 

40.00 

8.6686+05 

3 . 36 

324.4-j 

34.31 

?.O0 

2.790E+04 

1.50 

1 13.7.1 

20 . 00 

1 . 0006*02 

0.09 

100.00 

40.00 

7.6166+05 

7.24 

332.  15 

HO. 41 

20.  ?0 

2.7?0E+fl4 

1.50 

115.31 

20.00 

1 .0O0E+O2 

0 . 59 

100.00 

40.00 

S.245E+05 

7 . 09 

339.89 

106.. 32 

1 : .  d* 

2 . 7*06*04 

1.50 

1 13.31 

20.00 

1 . 0006-02 

0.09 

iOO. 00 

40.00 

5.037E+OS 

6.3" 

347.63 

101.17 

;2 .  ?0 

2.790E+04 

i  .50 

115.31 

20.00 

1.030E+C2 

0.0? 

100.00 

40 . 00 

4. 5 15E*55 

0.03 

355.33 

39.26 

13.0© 

.  2302 *05 

3.59 

115.95 

20.09 

1 . OO0E+02 

C .  05 

1 00 . 00 

40 . 00 

4. 7726+05 

6.20 

303.  12 

105.21 

14.30 

2. 1Z5£^0: 

4.  15 

1 15.45 

20.  CO 

1 .0eC6*O2 

0.09 

100.00 

40.  ca 

4.5566*05 

6 . 06 

370.S6 

88.02 

4 

0.00 

2. 1095+05 

4.22 

115.17 

20.00 

1.0006*02 

0 . 09 

100.00 

40.00 

3 . 3646+05 

5.21 

371.71 

Z3.56 

54.48 

• .  00 

2 .  2?  ■:£*►>  5 

4.70 

1 14.66 

20.OO 

1 . O0OE+O2 

0 . 09 

ICO. 00 

40 . 00 

2.940E+05 

4.87 

33 l . 20 

2?.  3  3 

-35.71 

2 . 00 

:.3i2£*05 

4 1  7  1 

114.51 

20.00 

1 .0006+02 

0 . 09 

1 00 . 00 

45 .  eo 

2. 6686+05 

4.64 

331.20 

3*.  67 

3 . 00 

2.2  296*05 

4.14 

115.35 

20.00 

1 . 0006*02 

0.0? 

100.00 

40 . 00 

2.3796*05 

4.34 

331.20 

■10 . 62 

4.00 

‘.05£c*05 

4.07 

116.  15 

20 .  CO 

1 . 0006-02 

0 . 29 

1 CO . 55 

45.00 

1 . 9576*55 

3.97 

•31.20 

-39. 6_ 

5 . 00 

.  3656+05 

3.33 

1 16.97 

20.0O 

1 .0006+02 

0 . 05 

100.00 

40 . 00 

1 . 803E+O5 

3.61 

33 1 . 20 

33.66 

6.  30 

.7 Owe *04 

1 . 50 

115.31 

2O.0O 

l  .5006*02 

5.05 

lvO.OO 

4C .  OP 

2.2126*05 

4.2-2 

304.85 

"5.63 

" .  00 

. 2326+04 

2.35 

120.50 

20.00 

1 .0006*02 

0.0? 

150.05 

4*3 . 00 

3.3416*05 

5.  1? 

246 . 33 

61.45 

?.oo 

.  51 2  E  *C5 

4.  00 

114. 10 

20.00 

1 .0006*02 

0  .  O* 

1 00 . 00 

40 . 50 

4.7842+05 

6.21 

247 .32 

62.22 

9.O0 

. 3 1 96+05 

3.33 

1 17.29 

20 . 00 

l  .0006+02 

0 . 09 

100.00 

40.09 

5. 51 46*05 

6.67 

255 . 57 

64.  33 

.639E+05 

7 . 67 

117.67 

20.00 

2.  3366*05 

4.74 

204. 1 l 

40.00 

5 . 307E+05 

6.54 

263.31 

67.  33 

; .  oo 

.5936+05 

3.53 

•  :  ~  '3Q 

20.00 

2.5086*05 

4 . 50 

137.05 

40.00 

6.  1326+05 

".03 

271.05 

6?.  33 

: .  O0 

.  *5035  *03 

7.33 

114.  15 

22.00 

2 . 6 19E+05 

4.60 

197.37 

40.00 

6.375E+85 

7.44 

273.30 

72.43 

3.30 

.3 12  £*05 

7.61 

114.  14 

20 .  CO 

2.674E+05 

4.£4 

192.5? 

45.00 

7.5876+95 

7.32 

256 . 54 

74.32 

4.00 

3 . 62 

114.05 

20.00 

2.6726*05 

4.64 

194.65 

40.09 

3. 166E+95 

8.  11 

234.29 

7" .  04 

3.00 

. 741 £*05 

3.64 

113.38 

20 .  OO 

2.6 156*05 

4.59 

199.66 

40 . 00 

3. JO  16 +05 

3.23 

302 . 02 

79.05  1 

*5 . 00 

. "OSt+CS 

3.71 

113.65 

20 . 00 

2 . 4506*05 

4.47 

208 . 5 1 

40.00 

3.2776*05 

8.  17 

303.7" 

30.23 

" .  00 

.861 £+05 

3.37 

113.53 

20.00 

8.6426*04 

2 . 64 

223.00 

40.00 

7.9596*05 

8.03 

317.51 

•  -  •  -  *  ii 

3.00 

.353E+05 

3.37 

1 -4.94 

20.00 

1 . 0006*82 

0.09 

100.00 

40.00 

7.5526*05 

7 .32 

325.25 

■33  •  27  | 

9 .  0© 

4.2925*04 

1.3" 

125.53 

22.00 

1 . 0006*02 

5 . 03 

100.00 

40 . 00 

3.2636+95 

3.16 

333.00 

33.17  2 

20 .  CO 

.  7502*04 

1.50 

115.31 

20 . 00 

i . 0006*02 

0 . 09 

100/00 

40.00 

8.2496+05 

•3.15 

340 . "4 

11 -.54  jz 

' 

: .  30 

. 05 1 c *05 

3 . 65 

116.73 

20 . 00 

1  .0006*02 

0.09 

105.55 

40 . 00 

b. "5 16+95 

" .  3c 

348.48 

114.02 

21.00 

. 5723*05 

3 . 23 

l  *6.23 

20.O0 

1 . 0006*02 

0.09 

100.00 

40 . 05 

5. 162E+05 

6.45 

356v 2 7 

36 . 38  9 

'3.00 

2. 0*55= t  05 

4.03 

1 15 . 66 

20.0O 

1 . 0006+02 

0.O9 

100.00 

40 .  C0 

4.066E+O5 

5.73 

363.97 

3 1  .  bb 

24.00 

2 . 2096*05 

4.22 

115.17 

20 . 00 

l  .0006*92 

5 . 09 

100.00 

40.05 

3. 3646+05 

5.21 

371.7: 

SS . 56  2 

5 

0 . 00 

. 60C£*05 

3 . 32 

116.33 

20.0O 

1 . 5006*02 

5.0? 

100.00 

40.00 

2.7646+05 

4.74 

379.53 

3?.??  4 

59.58 

1 . 00 

.751 £*05 

3 . 30 

l 16. 76 

20 . 00 

1 . 000E+O2 

0 . 09 

100.00 

45 . 00 

2 . 536E+05 

4.52 

37?.'4 

•39 .  *5 

-?0 . 14 

2 . 00 

.75; £+05 

3.74 

1  16.43 

20 . 00 

1 . 00OE+02 

0.09 

J00.00 

48 . 00 

2.2 106+05 

4.22 

379 . 74 

3?.7?  “* 

T  .  0C 

.  723E+05 

3.74 

1 16.57 

20 . 00 

1 . G00E*O2 

0.59 

100.00 

48 . 50 

1 .8286+05 

3 .  c3 

379.74 

S9.7S 

4.00 

- 707E *05 

3 .  "*! 

116. 55 

20.  JO 

1  .OOOE+OZ 

0.0? 

100.50 

40.00 

1. 5356+05 

3.5" 

373.74 

33.7?  ^ 

5 . 00 

.696E+05 

3.70 

117.35 

20 . 00 

1  . 0006*02 

O.0? 

100.00 

40.90 

1.6786+05 

3.63 

313.73 

73.36  | 

6 .  00 

. 244£*05 

4.25 

I 16.03 

20 .  CO 

1 .0OOE+O2 

0.09 

100.05 

40.00 

2 .  11 JE+05 

4.  13 

253.47 

63.51  9 

" .  00 

..  1646*05 

4.  13 

1 14.60 

20.00 

1 .0006*02 

0.0? 

100.00 

40 . 00 

2.9436+05 

4.  7? 

249 . 32 

62.33  | 

3.00 

2 . 0596+05 

4.07 

1 17.22 

20.00 

1.0006*02 

0.09 

100.00 

40.00 

3.349E+05 

5.57 

855.6? 

6J.33  } 

9 .  00 

.*  I* £*05 

3.23 

113. 13 

20.00 

1 . 00OE+O2 

0 . 09 

100.55 

40 . 00 

5. 2 1  “6+05 

6.49 

263 . 43 

67.31  ’ 

0 . 00 

•  £b  1E+05 

7 . 66 

113.50 

20 . 00 

2.3316-05 

4.  33 

203.31 

40.00 

6. 3316+05 

7 .  U-’ 

271.1? 

09.7?  , 

! .  00 

.6076*05 

7.60 

1 1 3.76 

20 . 00 

2. 4546+55 

4.  45 

1 ?? . 30 

40.00 

9.0376*05 

9.05 

278.92 

72 . 24 

2 . 00 

.  6086+05 

3 .  £0 

113.31 

20.00 

2.523H+05 

4.51 

197.78 

40 . 00 

3.7296*05 

8.3? 

286 . 66 

"4.01 

3 . 00 

. 7 22 £*05 

3 .  >2 

: I?.-? 

20.00 

2.5376+05 

4_  c  -> 

198.33 

40.00 

9 . 2636+05 

2.54 

294.40 

76 . 93 

4.00 

.554? -05 

3.65 

113.70 

20 . 00 

2 . 495E+05 

4.  4? 

202 . 32 

40.00 

9.0406*05 

0  _  C  ^ 

302.15 

79.67 

5 . 00 

. "273*05 

3.  "3 

113.56 

20 . 00 

2 . 3926+05 

4.33 

210.  19 

40.00 

9.3506*05 

3.31 

309. S3 

30.0? 

1-5.00 

. i 575*05 

7.  ?7 

113. 43 

20 .  C0 

3.7096+04 

2.65 

222 . 23 

40.00 

9.349E+05 

8.91 

317 .63 

32.  rc 

“ .  00 

. 508E+05 

7.62 

1 14. “2 

20 . 00 

l  .O00E+O2 

0 .  C9 

150.00 

40 . 00 

9.217E+05 

3.02 

325.33 

33.53 

3.00 

.3 l 2E+05 

3.22 

1 16.O0 

20 . 00 

1 .0006+02 

0.09 

IO0.OC 

40 . 00 

7.74SE+05 

7.90 

332.  12 

84.  5  ? 

?.O0 

.  77  ',£+05 

7 . 32 

116.5? 

20.00 

1 .0006*02 

0.03 

100.50 

40.00 

6.056E+O5 

6.?? 

340 . IS 

106.05 

20.00 

.535E+05 

3.62 

1 16.37 

20 . 00 

l . 00OE+O2 

0.55 

100.05 

40 . 05 

J. 87  IE *©5 

6.23 

749. 61 

86.27 

21.00 

.  P&5CE+05 

7.65 

116.67 

20.C0 

1.9006+02 

0 . 59 

1 00 . 00 

40.00 

4.070E+05 

5.73 

356.35 

37.  r 

.2 . 00 

. "3?E*05 

3.74 

1 18.46 

2O.O0 

1 .0O0E+02 

0 . 3? 

100.00 

40.00 

3.533E+05 

5.34 

364.0? 

33.06 

.3.00 

. 794E+05 

3.20 

116.35 

20.00 

1  .uOOE+02 

0.3? 

100.00 

40.00 

3. 1006 +<35 

5.O0 

2’1 .33 

88.3? 

'4.00 

. 3OOE+05 

7.32 

1 16. 33 

29.00 

1.0006+92 

0.03 

100.00 

40.00 

2.7S4E+05 

4.  74 

•79. 58 

3?.?? 

li 


45 


0 .  00 

1 . 7 1 -5*05 

113.04 

Table  15. 

20.00  1 . 0005+02 

(cont.) 

0.09  1OO.O0 

2. 38 75 *05 

373.49 

90.26 

.  ■’»? 

: . ? : ?7*05 

r .  50 

115. 04 

20.00 

I . 0C0E+O2 

€.00 

1 00 . GO 

40.00 

2.  S 3-35 *05 

4. 8: 

731.7- 

90.61 

-  .  JO 

1 .5  l?E*07 

3 . 50 

113.04 

20.00 

1 .0005*02 

0.O9 

100.00 

40 . 00 

2.3135+05 

4.  "2 

76 1 . 77 

90.51 

3 .  OO 

1 .  :* !  ?e  +0S 

3. 50 

1 17.04 

20.20 

1.CO0E+O2 

O.09 

100.00 

.m.co 

1 . S49t*G5 

3.96 

2S 1 . 77 

?O.70 

4.00 

1.5195+05 

2.30 

117. 04 

20 . 00 

l .  O00E+O2  • 

O.09 

IPO. 00 

40. C0 

l . 7266+05 

3 . 73 

TP  1.77 

90 . 5 : 

5 .  00 

-..519E+05 

3.30 

113.04 

20.00 

1 .OOOE*02 

G.09 

100.00 

40 . 00 

1.7935+05 

3. 7T 

"31.52 

i  1 . 6  3 

4 . 00 

1 .  c34£*t’5 

7  •  30 

122.88 

:o .  oo 

1 . 0005*07 

0.09 

100.00 

40.20 

2 . 050E+05 

4.09 

2-r, 74 

63.99 

7 . 00 

: .  22 '=>05 

4.00 

121 .98 

20.00 

1  . 0005*02 

0.09 

10O.C0 

40 .  OP 

2 . 436 E +05 

4.43 

253.53 

63.3? 

= .  00 

\ .  3*405*05 

t  _ 

114.29 

20 . 03 

I . JOOE+02 

0 . 09 

100.00 

40 .  GO 

3.0 122+05 

4.93 

254.60 

64. 50 

? .  JO 

*os 

T, 

1 1  •  .  . 3 

20.00 

1 . 0005*02 

0.73 

>•0.00 

40.00 

3.5775+03 

5.37 

262.35 

86 .  ?0 

0.00 

I .  “64e>05 

3. 77 

l 17. 25 

20.00 

1 . 0005*02 

r .  .:<? 

1 00 . 00 

40.00 

4. 1725+05 

5.50 

2 "0 . 09 

69. 71 

i .  :o 

1.6 7 16+05 

3.47 

1 3 . 7  7 

20 . 00 

2.2945*05 

4.:' 

208.94 

40 .  OO 

4.5126+05 

6.30 

277 .  ;  7 

r: 

: .  jo 

:  .5335*05 

3 .  £3 

113.94 

2O.0O 

2 . 3495*05 

4.  1? 

“07  is 

40  .  GO 

5 . S935+05 

6.35 

-ce,  _  ep 

"4. 

3  .  vO 

.  5  44c  *05 

3.54 

! ;3 . r5 

20.D0 

2.3375+05 

4. 39 

207.2? 

40.00 

9.3236+9-3 

7  . 

293.32 

76.2? 

4,  00 

1 .57  4c *03 

37T7 

1 1 2 . 92 

20 .  OO 

2 . 798c  *03 

4.  17 

208 .95 

40.00 

7.280E+05 

7766 

30 1 . G6 

7  i .  jo 

=  _ 

'•3-25*05 

3.  ’4 

1 13.34 

20.00 

2.2345*05 

4. 70 

212.3? 

•10.00 

1636*05 

7.60 

70S.. 38 

30.25 

i .  0  0 

’. .  3405*05 

7.85 

1 1 3 . 49 

20.09 

7.9:75*04 

2.45 

27 1 . 30 

40  .  OP 

9.7275*05 

7.37 

319.35 

?  t  .  ?4 

“ .  JO 

1 . 7? 16*05 

5.30 

1 14.42 

20.00 

i . 0005*0: 

P.03 

100.00 

40 . 00 

6.  2  996+03 

7.07 

324.2? 

•33.  40 

.  00 

1 . '045*03 

3 .  ~  1 

: 15.03 

20.00 

1 .0005*02 

0.03 

100.00 

40 . 00 

5.5236*0? 

8.7: 

332.07 

34.53 

0.00 

: . 5 f 0c *05 

37T5 

:  ;6.  75 

20.00 

1.0005*02 

0.09 

J 00.00 

40 . 00 

4.9845+05 

9.33 

339.76 

25.76 

0 . 00 

2.1375*03 

4.20 

i 12.25 

20.00 

1.000E+02 

0.0? 

100.00 

40 . 00 

4.4535+85 

3.29 

347.52 

36 . 45 

: .  oo 

2. J 375 *03 

4.20 

112.25 

20.00 

l .  OOOE+02 

0 . 0? 

:  00.00 

40.20 

4.0476+05 

3.7: 

755.29 

67.70 

2 .  ?0 

1.5135+05 

3.50 

113.04 

2O.00 

1.0005*02 

0  .  O? 

100-.  GO 

40 . 00 

3.6235*05 

5.40 

363.01 

33.16 

3 . 00 

1.5135*05 

7.50 

!  13.04 

20.00 

1.0005*02 

0.09 

100.00 

40.00 

3.2066+05 

5.06 

370.-5 

69,27 

4.00 

1.3135*05 

2.30 

113.04 

20.00 

1.0005*02 

0 .09 

•00.00 

40.7O 

2 . 6376+05 

4.22 

373.49 

90 . 26 

0.00 

1 .3135*05 

2.50 

l :  3 . 04 

20 .  OO 

: . 0005*02 

0 . 0? 

100. CO 

40 .  co 

2. ‘646*75 

4.72 

383. 15 

?)  .2? 

1 . 00 

1 .3135*03 

2.50 

117.04 

20 . 00 

1.0005*02 

0  .  O'3 

170.00 

40 . 00 

2 . 4536*05 

4.45 

737.15 

•  91.  16 

: .  oo 

;.  5  1.-5  *05 

3 . 50 

113.04 

20.00 

1.0005*02 

O.09 

1O0.O0 

40 . 00 

2. 1315+05 

4.  14 

- 

31  .  12 

?.oe 

1 . 515E+05 

7.50 

113. 04 

20.00 

1 .0005*02 

O.09 

1 00 . 00 

40.00 

1 . 3425+05 

3.96 

323 . 15 

?  •  .  i ; 

4.00 

» .5135*05 

3.50 

113.04 

20 .  OO 

1 .0OOE *02 

0.09 

ioo.ee 

40.  ec 

1.9S1E+05 

4.00 

742 . 32 

34.07 

•5 . 00 

1 . 5 1?E*05 

3,50 

113.04 

20 .  OP 

1 . 0005*02 

o 

100.00 

40.70 

2. 1396+05 

4.20 

TOO . 06 

75 . 44 

>3.00 

1. 3335*05 

7  .  =,c 

123.52 

20.00 

1 .0005*02 

0.03 

100.  *10 

4G .  20 

3 . 49  7  E +05 

4.46 

157 . 79 

65.1G 

■’.00 

2 .0 16E*05 

4.07 

12*3.88 

20.00 

1  .  OOOE+02 

0 . 0? 

:  oo . oo 

40 .  CO 

7 .  7?I5+<*5 

4.  74 

257.23 

65. i: 

2.00 

2.0535*05 

4.  10 

118.81 

20 .  OO 

1 .0005*02 

0.09 

100.00 

40 . 00 

3.  1 79E*05 

5.06 

259.27 

65.97 

3.  CO 

2.  133E+05 

4.  15 

1 16.33 

20.  GO 

l  .0005+02 

0.09 

100.00 

40.00 

3.61SE+05 

5.40 

267.01 

63.72 

0.00 

2.  180E*O5 

4.  1? 

1 15.34 

20 . 00 

1 . 3305+08 

“  .  75 

215.95 

40 . 00 

4 . 0995+05 

5.75 

274.  ,'5 

‘0.67 

: .  oo 

2.2 395 *03 

4.22 

1 14.30 

20.00 

2.  1 505+05 

4.  16 

211.7? 

40 .  oo 

4. 8 135+05 

6.10 

222.4? 

72.9? 

;  .oo 

2 . 2275*05 

4.22 

l i 4. o J 

20.00 

2.2095*03 

4.22 

210.33 

40  .  OO 

5. 235 E +05 

8.39 

290.24 

75.22 

3.00 

2.2225+05 

4.23 

1 14.34 

20 .  GO 

2 . 22  IE *05 

4.23 

211.42 

40.00 

5 . 3  12E+05 

6.34 

297.98 

‘7.34 

4.00 

2. 2045  *-05 

4.21 

1 15.53 

2O.O0 

2. 1855*05 

4.20 

215.  17 

40 . 00 

5.3366+05 

6 . 56 

305.72 

79.30 

5 .  oo 

2.  1705+05 

4 .  \i 

1  1  to  -  76 

20 . 00 

5 . 9*25+04 

2.14 

222 .02 

40 . 20 

5 . 233E+05 

6.50 

313.47 

81.09 

3 .  OO 

2. 1135*05 

a .  ;  ; 

113.5* 

20 . 00 

l . 0005*02 

0 .  '? 

100.00 

40.00 

3.063E+O5 

9.3? 

321.1! 

32 .  '0 

r .  oo 

2.O19E+05 

4.07 

120.39 

2  , .  OP 

1 . 3005*02 

.7 .  >09 

ICO. OP 

48 . 70 

4.9225*05 

6.23 

316.95 

84.  14 

5  .  ?0 

1 .34t£*05 

7. 38 

123.28 

20.00 

1. OOOE+02 

0 . 0° 

100.00 

40 . 30 

4. 560E*05 

6 . 06 

336 .8? 

55.23 

? .  00 

1.5135*05 

7.50 

113.04 

20.00 

i .  OOCE**?: 

0.09 

100.00 

40.20 

4.3565+95 

5.93 

244.44 

=e> . 

:o.oo 

2. 1975+05 

4.20 

i  1 2.25 

20 . 00 

1 . OOOE+02 

0.09 

1  CP . CO 

4C.00 

4.  13?E*03 

5.73 

352.  16 

37.22 

1 . 00 

2.  1375*05 

4.20 

U2. 25 

20.00 

1.0005*02 

0.09 

100.00 

40.00 

2 . 7926*05 

5.57 

759 . 92 

33.  10 

.00 

1.5 135+05 

3.50 

113.04 

20 .  OO 

1 . OOOE+02 

e.G? 

; oo . oo 

40.00 

3 .2‘3E+05 

r  .21 

io7 .  b.” 

33.  14 

:3.co 

1 .5135+05 

7.30 

: : 7. jd 

20.00 

I .0005+02 

0.09 

100.00 

40.00 

3.0306+05 

4.?J 

375.4 1 

90.13 

’4.00 

1 .5135*05 

7.50 

113. 04 

20.00 

•..•.o:e*?: 

-r . 

: : ?. go 

4.? .  ;o 

2. “946*0? 

4.72 

7:7.  1? 

?1 .23 

46 


3.4 


BURST  AND  CLOUD  OUTPUT 


I 


There  are  several  different  ways  to  obtain  printed  output  of 
cloud  (fireball,  debris  cloud  or  plume)  parameters.  Tables  16  through 
19  show  four  different  types  of  cloud  output:  1)  cloud  creation  and  split 
output  controlled  by  the  $BURSTS  print  flag  (Table  16);  2)  burst  summary 
output  provided  automatically  (Table  17);  3)  detailed  cloud  list  output 
controlled  by  the  $CLDOUT  output  flag  (Table  18);  and  4)  cloud  list  sum¬ 
mary  output  controlled  by  the  $CLDOUT  keyword  (Table  19). 

The  detailed  cloud  output  shown  in  Tables  16  and  18  share  a 
common  format.  The  labels  for  this  format  are: 


CLOUD 

- 

cloud  number  (assigned  internally  by  the  code) 

BURST 

- 

burst  number  (assigned  internally  by  the  code) 

TYPE 

type  of  cloud  (mixed  (ion/neutral)  fireball,  ionized 
fireball,  neutral  fireball,  late  plume  or  subsiding 

debris) 

TIME 

- 

time  that  data  is  valid  (GMT  in  modified  DHMS  format) 

LAT 

- 

latitude  of  cloud  center  (deg  N) 

LON 

- 

longitude  of  cloud  center  (deg  E) 

ALT 

- 

altitude  of  cloud  center  (km) 

RMAJ 

- 

major  horizontal  radius  (km) 

RMIN 

- 

minor  horizontal  radius  (km) 

RVERT 

- 

vertical  radius  (km) 

WXRAY 

- 

X-ray  yield  in  cloud  (MT) 

WFISS 

- 

Fission  yield  in  cloud  (MT) 

MASS 

- 

mass  of  cloud  (kg) 

VN 

- 

north  velocity  (m/s)  of  center  of  mass 

VE 

- 

east  velocity  (m/s)  of  center  of  mass 

VUP 

- 

upward  velocity  (m/s)  of  center  of  mass 

VXMAJ 

- 

major  axis  expansion  velocity  (m/s) 
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Table  16.  Cloud  creation  and  split  outout. 
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Table  17.  Burst  summary  output. 

■«*  BURST  SUMMARY  ** 


EL‘PST 

TIME 
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LATITUDE 
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LONGITUDE 
•:  r-EG'- 

ALTITUDE 
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.  MT ) 
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0.75 

2 

0: 

18:32:00 

40.00 

-30.00 

1. 00 

20.00 

10.00 

15.00 

Table  18.  Detailed  cloud  list  output. 
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Table  19.  Summary  cloud  list  output. 
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VXMIN  -  minor  axis  expansion  velocity  (m/s) 

VXUP  -  upward  expansion  velocity  (m/s) 

Af.IGN  -  fraction  of  ions  aligned  with  geomagnetic  field 

FRACV  -  fraction  of  debris  that  is  ionized 

MASSC  -  mass  in  conjugate  region  (kg)  (for  plumes  only) 

DTBURST  -  time  since  burst  (GMT  in  modified  DHMS  format) 

TBURST  -  time  of  burst  (GMT  in  modified  DHMS  format) 

TIMFLS  -  X-ray  emission  time  (GMT  in  modified  DHMS  format) 

LATFLS  -  X-ray  emission  latitude  (deg  N) 

LONFLS  -  X-ray  emission  longitude  (deg  E) 

ALTFLS  -  X-ray  emission  altitude  (km) 

DENFLS  -  atmospheric  density  (g/cc)  at  X-ray  emission  point 

DEC  -  declination  angle  (deg)  of  major  horizontal  axis 

The  labels  for  the  burst  and  cloud  list  summary  outputs  shown  in  Tables  17 
and  19  are  very  similar  to  the  above  labels  and  need  no  further  explanation. 

The  development  of  the  clouds  created  by  the  high  altitude  burst 
(burst  1  in  the  test  problem)  can  be  seen  in  Table  16.  The  burst  initially 
creates  a  mixed  (ion/neutral)  fireball  (cloud  1).  Approximately  2  minutes 
later,  after  rising  and  expanding  appreciably,  this  cloud  splits  into  an 
ionized  fireball  (cloud  3)  and  a  neutral  fireball  (cloud  4).  The  ionized 
fireball  in  turn  splits  into  a  late  plume  (cloud  5)  and  a  subsidence  cloud 
(cloud  6).  This  second  split  happens  over  13  minutes  after  the  burst,  well 
after  the  cloud  has  reached  apogee  and  started  falling.  Without  splitting 
or  changing  number,  the  neutral  fireball  (cloud  4)  turns  into  a  subsidence 
cloud  as  can  be  seen  in  Tables  18  and  19.  This  illustrates  that  a  single 
high  altitude  burst  can  create  many  clouds.  In  this  problem  the  high  alti¬ 
tude  burst  gave  rise  to  5  different  clouds  during  its  evolution  and  ended 
up  forming  3  different  clouds  (two  subsiding  debris  clouds  and  a  late  plume) 
at  late  times.  This  behavior  is  typical  of  most  high  altitude  bursts. 
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Low  altitude  bursts,  on  the  other  hand,  usually  display  less 
complicated  life  cycles.  The  typical  low  altitude  burst  creates  a  single 
neutral  fireball  which  rises  to  apogee  and  then  turns  into  a  subsidence 
cloud  on  the  way  back  down.  This  behavior  of  clouds  formed  by  low  altitude 
bursts  is  demonstrated  by  the  second  burst  in  the  test  problem.  The  sub¬ 
siding  debris  cloud  created  by  burst  number  2  (cloud  2)  is  seen  in  Tables  18 
and  19. 

3.5  HF  PROPAGATION  OUTPUT 

Like  most  HFNET  output  the  output  from  HF  propagation  calculations 
is  available  at  two  levels  of  detail.  The  more  detailed  output  is  generated 
during  the  execution  of  the  simulation  proper  and  is  optionally  produced 
depending  on  the  setting  of  the  $HFCALC  print  flag.  An  example  of  detailed 
simulation-time  Hr  propagation  output  is  shown  in  Table  20.  The  less  de¬ 
tailed  propagation  summary  output  is  produced  by  the  post -processor  program 
after  the  simulation  has  completed  execution.  An  example  of  post-processed 
propagation  summary  output  is  shown  in  Table  21. 

Table  20  shows  the  detailed  HF  propagation  output  for  the  Wash¬ 
ington  to  Cambridge  Bay  link  at  20:00  GMT.  (Detailed  output  for  the  other 
two  links  at  this  time  was  also  produced  during  this  run  but  is  not  shown.) 
The  first  few  lines  announce  the  event  as  an  "HF  CALCULATION  EVENT"  and  show 
the  values  of  the  $HFCALC  parameters  being  used.  Of  particular  interest 
are  the  frequency  of  -20  Megahertz  and  the  print  flag  of  1.  The  negative 
frequency  caused  the  MUF  propagation  calculation  to  be  suppressed  and  the 
print  flag  caused  the  detailed  propagation  output  seen  here  to  be  generated. 

Below  the  event  header  are  the  transmitter  and  receiver  node  names 
enclosed  in  a  box  of  asterisks.  Because  the  hops  parameter  on  the  $HFCALC 
data  line  was  ■'^ecified  as  zero,  the  code  computed  how  many  hops  were 
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Table  20.  Detailed  HF  propagation  output. 
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Table  20.  (cont.) 
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Table  21.  Summary  HF  propagation  output. 
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appropriate  for  this  link  using  the  great  circle  distance  from  the  trans¬ 
mitter  to  the  receiver.  In  this  case  the  G.C.  distance  is  3778  km  which 
translates  into  a  minimum  of  1  hop  and  a  maximum  of  4  hops.  Printouts 
from  the  1,  2,  3,  and  4  hop  propagation  calculations  are  shown  in  Table  20, 
separated  by  lines  of  equal  signs.  The  printouts  for  the  successful  1  and 
2  hop  modes  are  further  subdivided  into  three  subsections  labeled  "MODE 
GEOMETRY",  "PROPAGATION  LOSSES"  and  "SIGNAL/NOISE  CALCULATION". 

Within  the  mode  geometry  section  the  first  printout  has  the  fol¬ 
lowing  labels: 

F90  -  upper  decile  frequency  (MHz) 

F50  -  median  frequency  (MHz) 

F10  -  lower  decile  frequency  (MHz) 

FREQ  -  frequency  being  used  (MHz) 

PROB  -  probability  that  the  given  frequency  will  propagate 

These  quantities  give  an  indication  of  the  variability  one  can  expect  over 
time  in  the  ambient  ionosphere's  ability  to  support  the  given  mode.  The 
F90,  F50,  and  F10  frequencies  should  propagate  90%,  50%,  and  10%  of  the 
time  with  the  F50  frequency  being  the  N-hop  MUF  during  median  ionospheric 
conditions.  The  probability  that  the  given  frequency  will  propagate  is 
computed  from  the  three  decile  frequencies  assuming  a  chi-squared  distri¬ 
bution. 


In  some  cases  the  frequency  used  is  above  the  F50  frequency  (or 
MUF) .  Since  frequencies  above  the  MUF  have  a  reduced  chance  of  propagating 
during  median  ionospheric  conditions,  the  ionosphere  is  artifically  "en¬ 
hanced"  to  its  lower  decile  level  in  order  to  obtain  useful  information 
about  this  frequency  that  otherwise  would  be  lost.  Whenever  this  happens 
the  message  "**  USING  ENHANCED  IONOSPHERE  is  printed.  The  2,  3, 

and  4  hop  printouts  in  Table  20  all  display  this  message.  In  the  2-hop 
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case  the  frequency  (20  MHz)  is  slightly  above  the  2-hop  MUF  and  well  below 
the  2 -hop  lower  decile  frequency,  so  propagation  is  possible  with  the  iono¬ 
sphere  "enhanced"  to  the  10%  level.  However,  both  the  3  and  4  hop  lower 
decile  frequencies  are  below  20  MHz  so  the  3  and  4  hop  modes  fail,  even 
with  the  enhancement. 


Moving  along  to  the  second  printout  in  the  mode  geometry  section 
the  following  labels  are  seen: 


ELEVATION 

AZIMUTH 


TIME  DELAY 
RAY  PATH 


-  transmitted  and  received  elevation  angles  (deg) 

-  transmitted  and  received  azimuth  angles  (deg) 
(measured  relative  to  G.C.  azimuth  angle,  posi¬ 
tive  clockwise) 

-  group  time  delay  (ms) 

-  ray  path  distance  (km) 


Under  the  title  "INTERMEDIATE  RESULTS"  are  several  quantities 
which  show  the  intermediate  steps  in  the  mode  geometry  calculation.  The 
headings  are: 


POINT 

-  breakpoint  along  the  raypath  (point  1  is  the 
transmitter;  point  2  is  the  first  ionospheric 

virtual  reflection  point;  point  3  is  the  re¬ 
ceiver  for  1  hop  modes  or  the  first  ground 
reflection  point  for  multihop  modes;  etc.) 

EQUAL  SPACING 

latitude  and  longitude  (deg  N,  E)  of  equally 
spaced  points  along  the  great  circle  path  from 

transmitter  to  receiver 

SINGLE  HOP  V.H. 

-  virtual  reflection  heights  (km)  after  sequen¬ 
tially  accounting  for  the  effects  of  the  F2, 

FI  and  E  layers  at  each  hop 

EQUAL  ANGLES 

-  latitude  and  longitude  (deg  N,  E)  of  breakpoints 
after  adjusting  for  equal  angles  of  incidence 

and  reflection 
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ALONG  G.C. 


-  tilt  (deg),  displacement  (km)  and  new  virtual 
height  (km)  after  accounting  for  ionospheric 
gradients  parallel  to  the  great  circle  path 
ACROSS  G.C.  -  tilt  (deg)  and  displacement  (km)  due  to  ion¬ 

ospheric  gradients  perpendicular  to  the  great 
circle  path 

The  last  printout  in  the  mode  geometry  section  shows  the  final 
mode  geometry.  Under  the  title  "FINAL  AMBIENT  BREAKPOINTS"  are  the  headings 
"POINT",  "LAT",  "LON",  "V.H."  and  "LAYER"  which  need  no  further  explanation 
except  that  the  latitude  and  longitude  are  in  degrees  (north  and  east)  and 
the  virtual  height  in  kilometers. 


Immediately  after  the  "PROPAGATION  LOSSES..."  heading,  losses  in 
signal  power  are  shown  broken  down  into  four  categories  depending  on  the 
mechanism  which  caused  the  loss.  The  labels  are: 


FREE 

AMBIENT 

NUCLEAR 

OTHER 

TOTAL 


free  space  loss  (dB) 

losses  due  to  ambient  D-region  absorption  (dB) 
excess  nuclear  induced  absorption  in  the  D-region  (dB) 
losses  due  to  various  other  mechanisms:  ground  reflec¬ 
tion  and  scattering  losses,  deviative  losses,  etc.  (dB) 
total  losses  (dB)  (sum  of  the  four  component  losses) 


The  next  printout  in  Table  20  shows  details  from  the  ambient 
and  nuclear  D-region  absorption  calculations.  This  printout  is  complicated 
but  also  very  useful  in  the  interpretation  of  results.  The  labels  are  as 
follows : 

HOP  -  hop  number 

POINT  -  index  of  endpoint  for  this  pass  (point  2  indicates 

the  upward  pass  on  the  first  hop;  point  3  downward 
pass  on  the  first  hop;  etc.) 
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D-PASS  (DB) 


LAT 


LON 

SECANT 

%QGAM  CLD 

%QBET  CLD 

%QAMB 

NE  (Q) 

NE  (XR)  CLD 

ALT 


D-region  absorption  (dB)  during  this  pass  (the 
top  number  is  the  excess  nuclear  absorption; 
the  bottom  number  is  the  ambient  absorption) 
latitude  (deg  N)  of  ray  path  D-region  entry  and 
exit  points  for  this  pass  (the  top  number  is  for 
the  intersection  with  the  bottom  of  the  D-region 
and  the  bottom  number  is  for  the  intersection  with 
the  top  of  the  D-region) 

longitude  (deg  E)  of  ray  path  D-region  entry  and 
exit  points 

secant  factor  at  the  ray  path  D-region  entry  and 
exit  points 

percentage  of  total  Q  (ion  production  rate)  at 
altitude  ALT  due  to  the  largest  single  gamma 
source;  cloud  number  of  that  gamma  source 
percentage  of  total  Q  at  altitude  ALT  due  to  the 
largest  single  beta  source;  cloud  number  of  that 
beta  source 

percentage  of  total  Q  due  to  ambient  conditions  at 
altitude  ALT 

_3 

electron  density  (cm  )  due  to  Q  at  altitude  ALT 

-3 

electron  density  (cm  )  at  altitude  ALT  due  to  the 
largest  single  X-ray  source;  cloud  number  of  that 
X-ray  source 

altitude  (km)  of  the  largest  single  contribution 
to  the  total  absorption  on  this  pass 


As  an  example  of  the  utility  of  this  printout  in  interpreting 
results  refer  to  the  1-hop  and  2-hop  printouts  shown  in  Table  20.  The 
1-hop  mode  shows  an  ambient  D-region  absorption  loss  (AMBIENT)  of  3.82  dB 
and  an  excess  nuclear  D-region  absorption  loss  (NUCLEAR)  of  0.00  dB.  In 
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contrast  the  2-hop  mode  shows  an  ambient  loss  of  4.86  dB  and  an  excess 
nuclear  loss  of  15.55  dB.  Why  are  the  losses  so  much  higher  in  the  2-hop 
case? 


The  difference  in  the  ambient  absorption  is  easily  explained  by 
the  fact  that  the  2-hop  path  has  4  D-region  crossings  while  the  1-hop  path 
has  only  2.  Thus  the  total  absorption  is  greater  for  2-hops  than  for  1-hop 
even  though  each  pass  in  the  two  hop  case  has  less  dB  due  to  different 
secant  factors.  (1.09,  1.20,  1.25  and  1.32  dB  for  the  2-hop  path  versus 
1.79  and  2.03  for  the  1-hop  path.) 

The  difference  in  the  excess  nuclear  absorption  between  the  1-hop 
and  2-hop  modes  can  be  explained  by  two  observations.  The  first  is  that  the 
only  source  of  nuclear  absorption  shown  in  the  printout  is  cloud  2,  seen 
under  the  "%QGAM  CLD"  heading  for  both  modes.  (Cloud  2  is  a  subsiding  debris 
cloud  located  at  40.044°N,  -79.852°E  and  28.348  km,  according  to  the  print¬ 
out  in  Table  18.)  The  second  observation  is  that  the  geometry  for  a  1-hop 
mode  is  quite  different  from  that  of  a  2-hop  mode,  especially  in  the  loca¬ 
tion  of  the  D-region  crossings.  The  D-region  crossings  of  interest  in  this 
case  are  the  ones  closest  to  Washington.  As  shown  in  Table  20,  the  1-hop 
path  enters  the  D-region  at  43.25°N,  and  -79.01°E  (approximately  360  km 
from  the  cloud)  while  the  steeper  2-hop  path  enters  at  40.13°N  and  -77.64°E 
approximately  190  km  from  the  cloud).  Referring  back  to  Figure  1,  the 
cloud  is  approximately  in  the  same  location  as  the  burst,  marked  by  an 
asterisk,  and  the  1-hop  and  2-hop  D-region  entry  points  are  marked  with  a 
1  and  a  2,  respectively.  Now  it  is  easy  to  see  how  this  difference  in 
excess  nuclear  absorption  came  about,  and  also  how  mode  geometry  can  some¬ 
times  play  a  crucial  role  in  nuclear  effects  calculations. 
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Moving  on  to  the  signal  and  noise  calculation  printouts  (in  Table 
20  again)  we  find  the  following  labels  under  the  "NOISE  SOURCES  (DBW)"  title: 

THERMAL  -  receiver  thermal  noise  (dBW) 

GALACTIC  -  attenuated  galactic  noise  (dBW) 

MAN-MADE  -  man-made  (local)  noise  (dBW) 

ATMOSPHERIC  -  attenuated  atmospheric  (thunderstorm)  noise  (dBW) 
TOTAL  NOISE  -  total  noise  power  (dBW)  at  receiver 

These  noise  values  reflect  the  effective  transmitted  power  of  the  different 
noise  sources  as  well  as  propagation  losses  (for  galactic  and  atmospheric 
noise)  and  the  bandwidth  of  the  receiver.  To  the  right  of  the  noise  print¬ 
out  are  the  following  quantities: 

TRANS  POWER  -  transmitted  signal  power  (dBW) 

RECEIVED  POWER  -  received  signal  power  (dBW) 

(=  transmitted  power  -  propagation  losses) 

S/N  -  signal  to  noise  ratio  (dB) 

(=  received  power  -  total  noise) 

In  many  applications  the  signal  to  noise  ratio  is  the  "bottom  line",  but 
as  seen  above  many  of  the  quantities  provided  in  the  detailed  output  are 
needed  for  the  proper  interpretation  of  results. 

The  summary  HF  propagation  output  in  Table  21  contains  much  of 
the  same  information  as  the  detailed  output  in  Table  20.  One  thing  that 
should  be  noticed  however,  is  that  all  modes  calculated  for  this  link  at 
this  time  are  included  in  the  summary  output,  even  though  they  may  have 
been  calculated  during  separate  HF  calculation  events.  In  this  case  the 
modes  for  10  MHz  and  the  various  MUFs  are  shown  in  Table  21  in  addition 
to  the  20  MHz  modes.  The  summary  output  has  two  headings  that  require  some 
explanation: 


GO 


SYM  -  symbol  explaining  something  about  this  mode: 

*  =  computed  during  a  "multimode"  calculation 
M  =  this  is  the  N-hop  MUF 

E  =  computed  using  "enhanced"  ionosphere 
F  =  mode  failed  to  propagate 

MODE  -  number  of  hops  and  the  reflecting  layer  for  each  hop 
E  =  E-layer  reflection 

*  =  Fl-layer  reflection 
F  =  F2-layer  reflection 

3.6  PLUME  MODE  PROPAGATION  OUTPUT 

Output  from  the  plume  mode  calculation  is  also  available  at  two 
levels  of  detail.  Table  22  shows  the  detailed  plume  mode  output  generated 
during  the  plume  mode  calculation  at  20:00  GMT.  This  output  can  be  turned 
on  or  off  using  the  $PLCALC  print  flag.  Summary  plume  mode  output,  on  the 
other  hand,  is  produced  automatically  during  any  plume  mode  calculation 
event.  An  example  of  plume  mode  summary  output  is  shown  in  Table  23 
for  the  Santa  Barbara  to  Big  Fork  link  at  20:00  GMT. 

Referring  to  Table  22,  the  plume  mode  calculation  event  notice 
line  is  followed  by  printouts  of  the  plume  modes  for  each  link.  In  this 
case  two  plume  modes  were  found  for  the  Santa  Barbara  to  Big  Fork  link  and 
none  were  found  for  the  other  two  links.  The  following  quantities  appear 
on  the  first  line  of  output  for  each  successful  plume  mode: 

CLOUD  -  number  of  cloud  supporting  plume  mode 

TYPE  -  type  of  plume  mode;  three  types  are  defined: 

SIDE  -  side  mode  off  of  field  aligned  plasma 
BASE  -  forward  scatter  off  of  plume  base 
ISOT  -  isotropic  scatter  off  of  plume  base 

LAT  -  latitude  (deg  N)  of  reflection  point 
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Table  22.  Detailed  plume  mode  output. 
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Table  23.  Summary  plume  mode  output. 
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LON 

ALT 

V.H. 

FREQ 

CFREQ 

SIGMA 

DELAY 

ELT 

AZT 

ELR 

AZR 


longitude  (deg  E)  of  reflection  point 
actual  altitude  (km)  of  reflection  point 
virtual  height  (km)  of  reflection  point 
frequency  (MHz) 

critical  frequency  (MHz)  -  maximum  frequency  that  would 
reflect* 

2 

reflecting  cross  section  (km  )* 

group  time  delay  (ms) 

elevation  angle  (deg)  at  transmitter 

azimuth  angle  (deg)  at  transmitter  (measured  relative 

to  G.C.  path,  positive  clockwise) 

elevation  angle  (deg)  at  receiver 

azimuth  angle  (deg)  at  receiver  (measured  relative  to 
G.C.  path,  positive  clockwise) 


These  quantities  describe  the  basic  geometry  of  the  plume  mode. 


For  forward  scatter  base  modes  a  second  line  of  printout  is  also 
provided  which  describes  the  "footprint"  region: 


FPLAT 

FPLON 

FPRMIN 

FPRMAJ 

GRT 


latitude  (deg  N)  of  footprint  center 
longitude  (deg  E)  of  footprint  center 
minor  radius  (km)  of  footprint 
major  radius  (km)  of  footprint 

ground  range  (km)  from  footprint  center  to  receiver 


The  "footprint"  is  an  elliptical  region  on  the  surface  of  the  earth  illumi¬ 
nated  by  rays  propagating  from  the  transmitter  and  bouncing  off  of  the  base 
of  the  plume  (assumed  to  act  like  a  mirror).  If  the  receiver  is  in  or  very 
near  this  region  then  these  forward  scatter  plume  base  modes  are  predicted. 


* 

See  Reference  2  for  explanations  of  these  quantities. 
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The  remainder  of  the  plume  mode  detailed  output  consists  of  the 
propagation  loss  printout,,  the  detailed  D-region  absorption  printout  and  the 
noise,  signal  and  S/N  printout.  These  printouts  have  exactly  the  same  for¬ 
mat  and  meaning  as  their  counterparts  in  the  HF  propagation  output  explained 
above . 

As  can  be  seen  by  comparing  Tables  22  and  23,  the  plume  mode 
summary  output  contains  a  subset  of  what  the  detailed  output  contains, 
presented  in  a  more  concise  format.  Also  provided  in  the  summary  output 
are  link,  node  and  ambient  environment  parameters,  included  to  reduce  the 
amount  of  "page  flipping"  required  on  the  part  of  the  reader. 

3.7  SIMULATION  STOP  OUTPUT 

At  the  end  of  the  simulation  HFNET  reports  the  status  of  the  event 
list  and  dynamic  storage  in  the  same  format  as  at  the  beginning.  It  is 
interesting  to  compare  the  situation  at  the  end  of  the  test  problem  run 
(shown  in  Table  24)  with  that  at  the  start  (shown  in  Table  14) .  Note 
that  the  event  list  at  the  end  of  the  run  still  contains  events  to  be 
executed:  4  cloud  updates,  a  cloud  list  output  event,  2  HF  calculation  events 
and  a  plume  mode  calculation  event.  Also  note  that  dynamic  storage  has  a 
great  deal  more  structure  than  it  did  at  the  beginning  and  that  well  over 
30,000  calls  were  made  to  dynamic  storage  routines  during  the  run.  These 
observations  reveal  two  important  aspects  of  HFNET 's  simulation  structure. 

One  is  that  the  simulation  would  continue  executing  events  (mostly  cloud 
updates)  forever  if  it  weren't  for  the  stop  event.  The  second  is  that 
dynamic  storage  is  truly  the  heart  of  the  simulation  structure  from  which 
the  rest  of  the  program  obtains  its  sustenance  (events  and  data) . 

One  additional  printout  is  provided  at  the  end  of  a  run.  This 
is  the  event/module  runtime  statistics  table  seen  at  the  bottom  of  Table 
24.  In  addition  to  displaying  the  total  number  of  events  executed  and 
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the  amount  of  CPU  time  they  took  to  execute  (593  events  and  188.380  seconds 
in  this  case)  the  following  quantities  are  shown: 


EVENT 

TIMES 

MODULE 

CALLS 

CP  SEC 

MEAN 

MIN 

MAX 


event  number  (see  Appendix  A) 

number  of  times  this  event  was  executed 

module  number  (see  Appendix  A) 

number  of  times  this  module  was  called 

total  CPU  time  (sec)  used  by  this  event  or  module 

mean  CPU  time  (sec)  used  by  this  event  or  module 

minimum  CPU  time  (sec)  used  by  this  event  or  module 

maximum  CPU  time  (sec)  used  by  this  event  or  module 


This  printout  can  be  extremely  useful  in  understanding  the  workings  of  the 
simulation  and  also  in  locating  where  the  code  spends  most  of  its  time. 

In  this  case  the  majority  of  the  time  (131.600  seconds)  was  spent  executing 
HF  propagation  calculation  events  (event  4's).  As  can  be  seen  by  comparing 
this  table  with  the  tables  of  events  and  modules  in  Appendix  A,  almost  all 
of  the  defined  events  and  modules  in  the  program  were  exercised  during  this 
test  problem. 


3.8  DEBUG  OUTPUT 


Debug  output  is  optionally  produced  by  most  major  HFNET  modules 
and  subroutines.  This  type  of  output  is  only  produced  when  requested  by 
the  user  on  the  $DEBUG  data  line.  The  quantities,  units,  and  labels  are 
peculiar  to  the  routine  producing  the  output.  An  effort  has  been  made  to 
provide  ample  debug  print  statements  with  clear  labels  for  routines  where 
such  output  might  be  found  useful  in  the  location  and  extermination  of 
bugs.  However,  since  the  authors  of  the  code  suffer  from  human  limitations 
in  the  prediction  of  their  own  errors,  inadequacies  in  debug  output  can  be 
anticipated.  A  fairly  deep  understanding  of  the  workings  of  the  HFNET  code 
and  a  listing  of  the  FORTRAN  are  prerequisites  to  the  successful  interpre¬ 
tation  of  debug  output.  Naturally  the  creators  of  the  code  would  appre¬ 
ciate  hearing  about  any  bugs  found  and/or  fixed  by  users. 
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3.9 


TROUBLE  OUTPUT 


Many  subroutines  in  the  HFNET  code  make  tests  for  the  reasonable¬ 
ness,  consistency  and  existence  of  data  on  which  they  are  instructed  to 
operate.  Should  such  a  test  fail,  an  error  message  is  normally  printed  and 
the  program  aborted  by  a  call  to  subroutine  ABORT.  This  routine  then  prints 
a  message  telling  which  subroutine  requested  the  abortion  and  the  line  num¬ 
ber  in  that  routine  where  the  trouble  was  detected.  ABORT  will  also  attempt 
to  provide  additional  information  which  may  prove  useful  in  locating  the 
source  of  the  problem. 

In  addition  to  the  program  stopping  due  to  errors  that  are  dis¬ 
covered  and  reported  by  the  code  itself,  it  is  also  possible  that  the  pro¬ 
gram  execution  may  be  interrupted  due  to  errors  caught  by  the  hardware, 
the  operating  system  or  the  FORTRAN  run-time  system.  Problems  of  this  sort 
are  usually  the  result  of  programming  errors,  but  sometimes  are  due  to  bad 
input  data.  The  normal  course  of  events  when  this  happens  is  that  some 
sort  of  cryptic  error  message  is  supplied  by  the  system,  a  core  dump  and/or 
a  subroutine  call  traceback  is  provided,  and  the  job  is  terminated.  Of 
course,  in  theory,  if  the  computer  is  working  properly,  the  input  is 
reasonable,  and  no  bugs  reside  in  the  code,  no  such  problems  will  occur. 
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APPENDIX  A 


SIMULATION  STRUCTURE 


HFNET  is  a  forward -running  event -sequenced  simulation  that  employs 
a  flexible,  modular  logical  structure.  In  a  very  real  sense,  HFNET  consists 
of  several  independent  simulations  which  operate  under  the  control  of  an 
overall  simulation  ’’manager",  as  depicted  conceptually  in  Figure  A.l.  In 
this  chart,  the  ambient  environment  simulation  contains  descriptions  of  the 
location  and  characteristics  of  all  transmitters  and  receivers  as  well  as 
models  of  various  "ambient"  phenomena  such  as  RF  noise,  atmospheric  winds 
and  natural  ionospheric  conditions.  The  nuclear  environment  simulation 
includes  time-dependent  specification  of  the  location  and  properties  of 
debris  regions  formed  by  nuclear  bursts.  In  other  words,  the  ambient  and 
nuclear  environment  simulations  together  model  the  "real  world"  in  which 
the  propagation  simulation  must  operate.  Finally,  the  propagation  simulation 
computes  what  effects  this  modeled  environment  may  have  on  signals  propagat¬ 
ing  through  it. 

The  basic  HFNET  simulation  control  structure  is  illustrated  in 
Figure  A. 2.  As  shown,  the  simulation  consists  of  a  number  of  distinct  modules 
that  are  interconnected  by  the  simulation  manager.  A  module  is  a  subroutine 
(or  collection  of  subroutines)  which  is  a  self-contained  simulation  of  a 
specific  part  of  the  overall  simulation.  Modules  are  only  called  by  the 
simulation  manager,  never  by  any  other  module,  and  have  no  calling  arguments. 
All  necessary  inputs  are  obtained  from  the  event  data  list  and  from  appro¬ 
priate  dynamic  storage  lists  or  labeled  common  blocks.  The  modules  contained 
in  the  current  version  of  HFNET  are  listed  in  Table  A.l.  Modules  can  readily 
be  added  or  deleted  to  meet  the  requirements  of  each  application. 


Figure  A.l.  HFNET  conceptual  block  diagram. 
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Figure  A. 2.  HFNET  simulation  structure. 


Table  A.l.  HFNET  modules. 


Module  # 

Module  Description 

Event(s) 

1 

data  input 

1 

2 

simulation  start/stop  output 

1,  2 

3 

nuclear  burst  (cloud  creation) 

3 

4 

HF  propagation  calculation 

4 

5 

output  for  cloud  creation/update/split  events 

3,  6,  8 

6 

cloud  update 

6 

7 

compute  ionospheric  data  tables 

7 

8 

cloud  split 

8 

10 

HF/VHF  plume  mode  propagation 

10 

12 

end  of  run  housekeeping 

2 

13 

start  of  run  housekeeping 

13 

15 

read  wind  model  database 

15 

17 

output  for  ionospheric  data  tables 

7 

20 

output  cloud  data  for  plotting 

3,  6 

21 

cloud  list  output 

21 

22 

output  contour  plot  data 

22 

23 

cloud  list  input 

23 

(Modules  9,  11,  14,  16,  18,  19,  24,  and  25  are  available  for  future  use.) 


The  simulation  is  time-ordered  and  advances  by  means  of  "events" 
which  are  user  input  or  generated  internally.  A  HFNET  simulation  event  is 
defined  to  be  any  occurrence  that  requires  a  particular  sequence  of  calcula¬ 
tions  to  be  performed  at  a  specified  simulation  time.  For  example,  a  se¬ 
quence  of  calculations  which  constitute  the  event  "HF  propagation  calculation" 
might  include:  1)  calculation  of  the  maximum  useable  frequency  (MUF) ;  2) 
calculation  of  the  upper  and  lower  decile  frequencies;  3)  computation  of 
the  probability  of  mode  occurrence;  4)  determination  of  the  ray  path  geometry; 
5)  calculation  of  propagation  losses,  noise  power,  received  signal  strength 
and  signal-to-noise  ratio;  and  6)  printing  the  results.  For  structural 
purposes  these  calculations  are  considered  to  occur  simultaneously  even  though, 
strictly  speaking,  some  small  amount  of  "real"  time  does  pass  due  to  propa¬ 
gation  delays. 

Events  may  be  generated  in  any  time  order  by  any  module  and  are 
stored  in  the  event  list  which  is  processed  by  the  simulation  manager.  Each 
event  has  an  event  dataset  associated  with  it  which  contains  the  required 
information  to  process  the  event.  The  event  dataset  structure  is,  in  general, 
defined  differently  for  each  type  of  event.  However,  the  first  two  words 
of  every  event  dataset  specify  the  event  time  and  the  event  type.  Table  A. 2 
shows  the  event  types  defined  in  the  current  version  of  HFNET. 

HFNET  modules  and  HFNET  events  are  interconnected  by  means  of  an 
"event  transfer  list”.  The  event  transfer  list  is  operationally  a  part  of 
the  simulation  manager  and  specifies,  for  each  event,  the  sequence  in  which 
modules  are  to  be  called  to  execute  the  event.  This  module  transfer  list 
is  set  in  a  data  statement  and  can  readily  be  altered  to  redefine  an  event 
or  to  add  new  events  as  applications  dictate.  The  last  column  in  both 
Tables  A. 1  and  A. 2  show  these  event/module  interconnections  and  are,  in 
fact,  representations  of  the  module  transfer  list. 
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Table  A. 2.  HFNET  events. 


Event  # 

Event  Description 

Module ( 

1 

start  simulation/read  input 

1,  2 

2 

stop  simulation 

2,  12 

3 

nuclear  burst/cloud  creation 

3, 

5,  20 

4 

HF  propagation  calculation 

4 

6 

cloud  update 

6, 

5,  20 

7 

initialize  ionospheric  model 

7,  17 

8 

cloud  split 

8, 

5,  20 

10 

HF/VHF  plume  mode  propagation  calculation 

10 

13 

dynamic  storage  cleanup 

13 

15 

initialize  wind  model 

15 

21 

cloud  list  output 

21 

22 

contour  plots 

22 

23 

cloud  list  input 

23 

(Events  5,  9, 

11,  12,  14,  16-20,  24  and  25  are  available 

for  future  use 
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The  HFNET  simulation  structure  is  set  up  to  allow  individual  modules 
to  perform  calculations  on  a  time  scale  which  is  consistent  with  the  require¬ 
ments  of  the  objects  or  processes  being  simulated.  Interpolation  procedures 
are  used  to  provide  this  flexibility.  Data  computed  and  stored  by  a  module 
at  time  intervals  determined  by  it's  own  accuracy  requirements  can  then  be 
accessed  and  interpolated  by  other  modules  operating  at  higher  (or  lower) 
simulation  rates.  The  primary  example  of  this  in  HFNET  is  the  interaction 
between  the  nuclear  environment  simulation  and  the  propagation  simulation. 

The  user  can  specify  the  location,  time,  yield,  etc.,  of  any  number 
of  nuclear  bursts.  Once  the  burst  event  has  been  executed  and  the  nuclear 
debris  cloud  (fireball)  created,  the  model  itself  decides  the  rate  at  which 
the  cloud  data  need  to  be  updated.  Cloud  data  are  computed  and  stored  in 
the  cloud  list  at  two  simulation  times  which  always  bracket  the  current 
simulation  time.  Two  different  clouds,  or  the  same  cloud  at  different  points 
in  its  development,  are  generally  updated  at  different  rates.  For  example, 
low  altitude  bursts  usually  require  less  updating  than  high  altitude  ones 
and  cloud  phenomenology  at  late  times  can  normally  be  updated  at  a  slower 
rate  than  at  early  times. 

However,  the  propagation  simulation  does  not  know  (or  care)  about 
these  details  in  the  nuclear  environment  simulation.  What  the  propagation 
model  needs  to  know  is  the  location  and  parameters  of  all  existing  nuclear 
debris  regions  at  the  time  of  an  HF  skywave  or  plume  mode  calculation.  An 
interpolation  subroutine  (INTERP)  is  provided  to  supply  this  information. 

This  enables  the  propagation  model  to  work  independently  of  the  nuclear 
environment  model  and  greatly  reduces  the  impact  on  one  model  of  modifica¬ 
tions  to  the  other. 

Most  of  the  data  in  the  HFNET  simulation,  either  input  data  or 
internally  generated  data,  is  stored  in  a  flexible  dynamic  storage  allocation 
system.  This  permits  large  data  arrays  to  be  handled  efficiently,  and  avoids 
the  necessity  of  dimensioning  data  arrays  for  the  largest  number  of  allowed 
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quantities.  The  dynamic  storage  allocation  system  was  developed  at  MRC  and 
is  fully  documented  in  Reference  3.  It  provides  a  versatile  and  easy-to-use 
capability  for  data  storage  and  retrieval  of  dynamically  varying  data  lists 
and  list  structures.  The  dynamic  storage  data  lists  used  by  the  current 
version  of  HFNET  are  shown  in  Table  A. 3  below. 


Table  A. 3.  HFNET  dynamic  storage  data  lists. 


LIST  DATASET 

NAME  SIZE  DESCRIPTION 


EVNT 

NODE 

LINK 

IONO 

CLOD 

NOIS 

PLUM 

F77X 

WDML 

WDMH 

WDDI 


16  event  list 

20  node  list 

10  link  list 

1576  link  ionospheric  data  tables 

51  cloud  list 

50  noise  calculation  scratch  storage 

55  plume  mode  calculation  scratch  storage 

183  geomagnetic  data  temporary  storage 

9  median  wind  model  low  altitude  database 

33  median  wind  model  high  altitude  database 
87  diurnal  wind  model  data  base 
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APPENDIX  B 

ROUTINES  AND  COMMON  BLOCKS 


The  HFNET  computer  code  consists  of  a  main  program,  a  block  data 
routine,  195  subroutines  and  functions,  and  40  common  blocks.  The  following 
sections  contain  alphabetical  lists  of  these  program  components  along  with 
a  brief  explanation  of  each. 

B.l  LIST  OF  HFNFT  ROUTINES 

ROUTINE  DESCRIPTION 

ABORT  Provides  a  subroutine  traceback  and  stops  the  program 

whenever  a  fatal  error  is  detected  by  any  HFNET  routine. 

ABSORB  Computes  incremental  D-region  absorption  given  the 

frequency,  electron  density  and  collision  frequency. 

ACLRX  Clears  an  array.  Alternate  entry  points:  ACLRR  and 

ACLRI . 

ADDCLD  Adds  a  newly  created  cloud  of  any  type  to  the  cloud 

list . 

AMOVX  Moves  an  array.  Alternate  entry  points:  AMOVR  and 

AMOVI. 

ASETX  Sets  an  array  to  a  specified  constant  value.  Alternate 

entry  points:  ASETR  and  ASETI. 

ASWICH  Keeps  track  of  several  logical  switches  relating  to 

the  mass  entrainment  parameter  of  a  rising,  expanding 
fireball . 
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ATMOSU 


A  VWIND 

BDRAG 

BEND 

BLACK 

BLKDAT 

BOUNCL 

BSCATR 

CHI 

CHISQ 

CLIP 

CLIPM 

CM2  KM 
COLUMN 


Provides  mean  atmospheric  quantities  as  a  function  of 
altitude  in  an  undisturbed  atmosphere. 

Finds  the  average  wind  velocity  at  the  center  altitude 
of  a  debris  cloud  by  averaging  over  the  vertical  thick¬ 
ness  of  the  cloud. 

Decelerates  the  velocity  component  of  a  cloud  perpen¬ 
dicular  to  the  geomagnetic  field  line. 

Calculates  the  amount  of  deflection  a  ray  suffers 
while  passing  through  an  ionized  layer. 

Black  body  deposition  function. 

Block  data  routine  used  to  perform  compile-time  ini¬ 
tialization  of  HFNET  common  blocks. 

Computes  ground  reflection  and  scattering  losses  as 
a  function  of  takeoff  angle,  frequency,  and  surface 
type  (land  or  sea) . 

Geometry  routine  for  plume  side  scatter  modes. 

Computes  solar  zenith  angle  as  a  function  of  position, 
time  and  season. 

Calculates  cumulative  chi-square  probability  function. 
Used  to  compute  mode  probability  as  a  function  of  the 
MUF  and  upper  and  lower  decile  frequencies. 

1-D  clipping  function  (real). 

Generalized  modulus  function.  Finds  X  modulo  I  where 
I  =  [A,  B]  is  any  interval. 

Convert  centimeters  to  kilometers. 

Computes  integrated  column  density.  Used  in  calcula¬ 
tion  of  X-ray  flash  electron  density  and  gamma  ioni¬ 
zation  rate  (QGAM) . 
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COMB 

COMPUT 


CPATCH 


CPUSEC 

CR2D 

CREATF 

CROSPT 

CUBIC 

D2CR 

D2R 

DAMBNT 


DATAPT 


Used  to  compute  all  possible  modes  for  multipath 
propagation  calculation. 

Used  by  the  ITS  and  RADC  (polar)  models  to  compute  F2 
layer  critical  frequencies  from  the  ITS  annual  and 
solar  cycle  coefficients.  Alternate  entry:  COMPT. 

Computes  parameters  for  the  creation  of  a  pseudo¬ 
fireball  in  the  opposite  hemisphere  (the  conjugate 
region)  to  simulate  the  debris  energy  patch  in  that 
region. 

Returns  CPU  time  used.  This  routine  is  specific  to 
the  VAX/VMS  operating  system. 

Convert  co-radians  to  degrees. 

Computes  initial  conditions  for  a  mixed,  neutral  or 
ionized  fireball  given  burst  location  and  yield. 

Subroutine  used  to  find  D-region  crossing  points  and 
secant  factors. 

Solves  cubic  equations.  Used  by  the  RADC  (polar) 
model  in  fitting  electron  density  profiles. 

onvert  degrees  to  co-radians. 

Converts  degrees  to  radians. 

Determines  ambient  D-region  electron  density  as  a 
function  of  altitude,  position,  solar  zenith  angle 
and  E-layer  critical  frequency. 

The  major  interface  between  the  nuclear  environment 
section  and  the  propagation  section.  Computes 
electron  density  and  collision  frequency  at  D- 
region  crossing  points. 
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DCROSS 

DEBUG 

DELCLD 

DELIST 

DELSET 

DENS 

DHMS2H 

DHMS2S 

DISCAT 

DISPLC 

DISTNC 

DPROB 

ECOMB 

F.INT1 

ELCOL 


Locates  all  D-region  crossing  points  for  a  ray-path. 

Tells  whether  or  not  debug  print  has  been  requested 
from  a  given  routine. 

Deletes  a  specified  cloud  from  the  cloud  list. 

Dynamic  storage  routine  to  delete  an  entire  D.S.  list. 

Dynamic  storage  routine  to  delete  a  specified  D.S. 
dataset . 

A  function  used  in  RADC  model  of  altitude  of  maximum 
electron  density  in  the  F2  layer. 

Converts  from  real  DHMS  format  to  hollerith  DHMS  format. 

Converts  real  DHMS  format  to  time  in  seconds. 

Calculates  distances  in  the  special  geomagnetic  XYZ 
coordinate  system  used  in  the  plume  side  mode  geometry 
calculation. 

Computes  displacements  due  to  ionospheric  tilts  (both 
parallel  and  perpendicular  to  the  great  circle  path.) 

A  function  used  by  DISCAT. 

Computes  probability  of  mode  occurrence  using  the 
MUF  and  the  upper  and  lower  decile  frequencies. 

Determines  E-layer  parameters  from  solar  E  and  sporadic 
E.  Part  of  the  RADC  (polar)  ionospheric  model. 

Exponential  integral  of  the  first  kind.  Used  in 
D-region  electron  density  calculation. 

Computes  the  location  of  point  2  given  the  location 
of  point  1  and  the  azimuth  angle  and  great  circle 
distance  of  point  2  from  point  1. 


84 


ELLIPS 

EQUATR 

EXOT 

EXPND 

F2DIS 

FEX 

FGENRL 

FLASHC 

FOEFUN 

GEOPT 

GETSET 


Computes  new  ellipse  radii  for  an  ellipse  being  per¬ 
turbed  by  wind  gradients.  Part  of  the  late  time 
subsiding  debris  model. 

Subroutine  to  compute  the  great  circle  distance  to 
the  equator  from  any  point  on  the  earth  given  the 
location  of  the  point  and  an  azimuth  angle.  Part  of 
the  atmospheric  noise  model. 

A  function  called  by  SCALHT.  Used  in  determining  the 
F2  layer  altitude  for  the  RADC  (polar)  model . 

Computes  the  radial  expansion  of  a  mixed,  ionized, 
or  neutral  fireball  during  a  single  time  step. 

Computes  upper  and  lower  decile  frequencies  given  the 
MUF,  sunspot  number,  location  of  point  and  time.  This 
routine  uses  data  from  the  DECILE.BIN  data  file. 

Computes  a  function  used  in  determination  of  the  radial 
expansion  of  a  fireball. 

A  function  used  in  obtaining  auroral  region  E  layer 
critical  frequencies;  called  by  QCHAT. 

Computes  X-ray  flash  electron  density  given  initial 
value  of  electron  density  and  rate  constants. 

Computes  E  layer  critical  frequency  as  a  simple  analy¬ 
tic  function  of  sunspot  number  and  solar  zenith  angle. 

Converts  from  geomagnetic  (dipole  field)  to  geographic 
coordinates.  This  is  the  inverse  routine  of  subroutine 
MAGNET. 

Dynamic  storage  routine  to  retrieve  a  specified  D.S. 
dataset . 
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GTRANS 

HEAVIS 


HEIGHT 

HFABSB 

HFMODE 

HFNET 

HFPOST 

HFS2NR 

HZ2MHZ 

ICLIP 

ICLIPM 

INPERR 

INTERP 

INTWND 

IONFIT 


Computes  atmospheric  gamma  energy  transmission  fraction 
given  the  column  density.  Used  in  calculation  of 
gamma  ionization  rate  (QGAM) . 

Interpolates  ionospheric  parameters  (layer  critical 
frequency,  altitude  and  semithickness)  for  three  ion¬ 
ospheric  layers  (E,  FI  and  F2)  from  precomputed  iono¬ 
spheric  data  tables.  Used  by  the  propagation  models. 

Computes  virtual  height  of  rays  reflected  by  parabolic 
ionospheric  layers. 

Computes  absorption  and  other  propagation  losses  on 
an  HF  skywave  ray  path. 

Computes  ray  path  geometry  for  HF  skywave  modes. 

Main  program. 

Writes  propagation  data  to  post  processor  output  file. 

Subroutine  to  calculate  and  print  noise  power,  signal 
strength  and  signal -to-noise  ratio  for  an  HF  skywave 
ray  path. 

Converts  hertz  to  megahertz. 

1-D  clipping  function  (integer). 

Generalized  modulus  function  (integer).  Finds  X  modulo 
I  where  I  =  [A,  B]  is  any  interval. 

Reports  input  errors  found  by  the  input  module  (MODI) . 

Computes  interpolated  cloud  data  for  a  specific  cloud 
at  a  specific  time.  Interpolated  cloud  data  is  stored 
in  common  block  NOWDTA. 

Computes  wind  dataset  index  given  the  altitude. 

Computes  best  fit  parabolic  parameters  for  Aerospace 
ionospheric  model. 
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IONUP  Returns  ionospheric  data  at  any  of  seven  points  between 

transmitter  and  receiver  interpolated  to  a  specified 
local  time. 

ITSION  Computes  ionospheric  parameters  for  24  hours  above 

a  specified  point  using  ITS  ionospheric  model . 

ITSPAR  Calculates  ionospheric  data  tables  for  ITS  ionospheric 
model . 

KM2CM  Converts  kilometers  to  centimeters. 

LAYER  Routine  used  to  compute  atmospheric  tilts. 

LISDAT  Dynamic  storage  routine  to  initialize  constants  and 

variables  in  LISTOR  common. 

LOSPRP  Checks  for  line-of-sight  propagation  using  a  4/3  earth's 

radius  to  account  for  atmospheric  refraction. 

LSWICH  Keeps  track  of  the  mass  entrainment  parameter  of  a 

rising,  expanding  fireball  relative  to  several  logical 
switches. 

LTRAN1  Converts  cloud  radii  from  a  coordinate  system  aligned 

with  the  geomagnetic  field  lines  to  one  aligned  with 
the  local  vertical. 

LTRNV1  Transforms  cloud  radial  (expansion)  velocities  from 

a  coordinate  system  aligned  with  the  geomagnetic  field 
lines  to  one  aligned  with  the  local  vertical. 

LTRNV2  Transforms  cloud  radial  (expansion)  velocities  from 

a  coordinate  system  aligned  with  the  local  vertical 
to  one  aligned  with  the  geomagnetic  field  lines. 

MAGFIN  Computes  magnetic  field  components  using  normalized 
magnetic  field  coefficients.  Used  by  the  ITS  iono¬ 
spheric  model. 
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MAGNET 

MAPOUT 

MAPSET 

MEMLOC 

MHZ2HZ 
MODI 
MODIO 
MODI  2 

MOD13 

MOD15 

MODI  7 

MOD  2 

M0D20 

MOD21 

MOD22 

MOD22X 

MOD23 


Computes  magnetic  field  components  using  an  earth 
centered  dipole  field.  Used  by  most  models  in  the 
code  requiring  magnetic  field  data. 

Dynamic  storage  routine  to  print  D.S.  maps. 

Initializes  map  projection  data  in  MAPCOM  common. 

Dynamic  storage  routine  to  compute  relative  address 
in  D.S.  buffer  given  absolute  address.  Also  performs 
buffer  swapping  to  secondary  storage  if  necessary. 

Converts  megahertz  to  hertz. 

Input  module. 

Module  to  control  plume  mode  calculations. 

Module  to  write  "end  of  simulation"  data  to  post  pro¬ 
cessor  data  file. 

Module  to  delete  F77X  data  list  to  free  dynamic  stor¬ 
age  space. 

Module  to  initialize  wind  model  data  lists. 

Module  to  print  ionospheric  data  tables. 

Module  to  print  simulation  start/stop  output. 

Module  to  write  debris  model  binary  (plot)  output. 

Module  to  generate  cloud  list  print  and  binary  output. 

Module  to  write  data  file  for  contour  plots. 

Used  in  conjunction  with  MOD22  to  compute  data  for 
"special"  contour  plots. 

Cloud  list  input  module.  Reads  in  precomputed  cloud 
data. 
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MOD  3 


MOD4 

MOD5 

M0D6 

MOD7 

MOD8 

MUFS 

MULMOD 

MULMUF 

MULTPH 

NINT 

NOISE1 

NOISE2 


Cloud  creation  (burst)  module.  Controls  the  creation 
and  storing  of  initial  fireball  parameters.  Also  stores 
an  update  event  for  newly  created  cloud. 

Module  to  control  HF  calculation  events. 

Output  module  for  cloud  creation,  split  and  updates. 

Cloud  update  module.  Controls  the  updating  of  all 
five  types  of  clouds. 

Ionospheric  initialization  module.  Controls  the  com¬ 
putation  of  link  ionospheric  data  tables  using  the 
three  ionospheric  models. 

Cloud  split  module. 

Computes  the  N-hop  maximum  useable  frequency  (MUF) . 

Computes  propagation  mode  geometry  for  "multimode" 
calculation. 

Computes  frequencies  to  use  during  "multimode"  cal¬ 
culations. 

Adjusts  ray-path  geometry  for  differing  reflection 
heights  and  ionospheric  tilts  during  "multimode" 
calculation. 

Nearest  integer  function.  Supplied  for  systems  that 
lack  this  routine.  Not  used  in  VAX  version. 

Used  by  the  noise  model  to  precompute  quantities  used 
in  the  atmospheric  noise  calculation  for  a  specific 
receiver  at  a  specific  time.  This  information  is 
saved  in  the  NOIS  data  list  for  subsequent  use. 

Computes  thermal,  galactic,  man-made,  and  atmospheric 
noise  power  for  a  specified  receiver,  time  and  frequency. 
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NXEVNT 


Searches  the  event  list  for  the  event  with  the  smallest 
time. 


NXTCRD 

NXTSET 

OB 

OPENFILES 

PARABI 

PATCH 


PATHW 

PLABSB 

PLMODE 

PLS2NR 

POLAR 

POLPAR 

PRABSB 


Routine  to  skip  comment  lines  in  the  input  file. 

Dynamic  storage  routine  used  to  cycle  through  a  D.S. 
data  list. 

Logical  function  to  test  for  out  of  bounds  variables. 

Routine  to  open  data  files  for  readonly  access.  Allows 
sharing  of  these  input  files  in  the  VAX  version. 

Computes  ionospheric  parameters  for  the  Aerospace 
ionospheric  model . 

Moves  the  burst  point  down  the  magnetic  field  line  to 
a  point  where  the  energy  trapped  in  the  local  conjugate 
is  a  certain  proportion  of  the  total  yield.  Most  of  the 
energy  and  debris  will  be  deposited  near  this  point. 

Routine  used  by  the  propagation  routines  to  locate 
data  in  an  ionospheric  data  table.  Also  computes 
spatial  interpolation  weights. 

Computes  propagation  losses  for  a  plume  mode  ray  path. 

Calculates  plume  modes  for  a  single  link  and  cloud. 

Computes  noise,  signal  and  signal-to-noise  ratio  for 
a  plume  mode  ray  path. 

Function  used  by  the  Aerospace  ionospheric  model. 

Computes  ionospheric  parameters  for  the  RADC  (polar) 
ionospheric  model . 

Generates  detailed  printout  of  propagation  losses  on 
an  HF  skywave  ray  path. 
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PREPAR 


Initializes  for  the  RADC  (polar)  ionosphere.  Reads 
and  stores  data  from  the  ITS  yearly  data  file  (YEAR.BIN) 
and  the  geomagnetic  data  file  (GEOMAG.BIN). 

PRMODE  Prints  HF  skywave  mode  geometry. 

PRMULT  Prints  "multimode"  output. 

PRNCLD  Prints  cloud  parameters  for  a  specified  cloud  at  a 

specified  time. 

PRNEVN  Generates  a  printout  of  the  event  list. 

PR0F3  Computes  a  function  used  by  PROFIL. 

PROFIL  Computes  vertical  electron  density  profiles  for  the 

RADC  (polar)  model. 

PRPLUM  Prints  plume  mode  geometry,  losses,  signal  and  noise 

data.  Also  writes  plume  mode  data  to  the  postprocessor 
data  file. 

PSTATS  Prints  event/module  runtime  statistics  at  end  of  run. 

QAMBNT  Computes  ambient  D-region  ion  production  rate  (Q) . 

QCHAT  Computes  ionospheric  parameters  derived  from  the 

Chatanika  electron  density  profiles.  Part  of  the 
RADC  (polar)  model. 

QEPHEM  Computes  solar  declination  and  ascension. 

QF2FH  Calculates  F2  layer  critical  frequency  for  RADC  (polar) 

model . 

QFETCH  Fetches  geomagnetic  data  stored  in  the  F77X  data  list 

for  use  in  computing  corrected  geomagnetic  coordinates. 
See  ZCGMCS.  Part  of  the  RADC  (polar)  model. 

QF0EF1  Computes  E  and  FI  layer  critical  frequencies  for  the 

RADC  (polar)  model. 
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QHTF2 


QMAGFI 

QOLION 

QUADR 

QUASIC 

R2D 

RAB 

RATES 

RCLCLD 

REFLC 

REMTABS 

RINTRP 


Computes  the  height  of  the  F2  layer  for  the  RADC 
(polar)  model. 

Computes  magnetic  field  components  for  the  RADC  (polar) 
model.  Same  as  MAGFIN. 

Computes  ionospheric  parameters  using  the  RADC  (polar) 
ionospheric  model . 

Solves  quadratic  equations  using  the  quadratic  formula. 
Used  by  the  RADC  (polar)  model. 

Computes  the  electron  density  and  collision  frequency 
at  a  D-region  point  given  the  pressure,  air  density, 
ion-production  rate  (Q)  and  the  attachment/detachment 
rates . 

Converts  radians  to  degrees. 

Computes  the  ground  range  and  azimuth  angle  from  point 
1  to  point  2  given  the  locations  (latitude,  longitude) 
of  the  two  points. 

Computes  D-region  attachment  and  detachment  rates. 

Recalls  cloud  data  for  a  specified  cloud  from  the 
cloud  list. 

Computes  the  maximum  frequency  that  will  reflect  off 
a  plume/fireball  and  the  effective  reflecting  cross 
section . 

Removes  tab  characters  (replacing  them  with  the  appro¬ 
priate  number  of  blanks)  from  80  character  card  images. 
Used  in  the  VAX  version  of  HFNET  to  make  echoed  input 
lines  print  correctly. 

Two  dimensional  (4  point)  interpolation  routine  used 
by  the  wind  models  to  compute  wind  velocities  at  the 
location  of  the  cloud  center. 


RISE 

S2DHMS 

SCALHT 

SCAT 

SETKP 

SETLNK 

SETPAR 

I 

SIMMGR 

SMOOTH 

SPACE 

SPLIT1 


Calculates  the  buoyant  rise  and  entrainment  of  ambient 
air  by  a  rising,  expanding  fireball. 

Converts  time  in  seconds  to  DHMS  format. 

Subroutine  used  by  the  RADC  (polar)  ionospheric  model 
to  obtain  altitude  of  maximum  electron  density  in  F2 
1 ayer . 

Subroutine  used  by  the  plume  mode  geometry  model. 

Computes  K  (corrected  Kp  index)  as  a  function  of  Kp 
index,  corrected  geomagnetic  time  and  latitude.  Part 
of  the  RADC  (polar)  ionospheric  model. 

Initializes  for  a  specified  link  in  preparation  for 
an  HF  propagation  calculation. 

Subroutine  to  initialize  constants  and  read  in  data 
for  the  ITS  ionospheric  model.  Data  is  read  from 
the  ITS  yearly  and  monthly  data  files  (YEAR.BIN  and 
MONTH. BIN) . 

Simulation  manager.  Manages  the  execution  of  all 
events  in  the  simulation  and  collects  data  on  event/ 
module  activity. 

Function  with  smooth  derivatives. 

Smooth (X)  =  j  2X2  O^X  ^.5 

j  1-2 (1-X) 2  .5<  X  ^1 

Function  used  by  the  wind  interpolation  routine  (AVWIND) 
to  compute  the  vertical  spacing  between  data  points 
in  the  wind  database. 

Handles  the  splitting  of  a  mixed  fireball  into  a  neutral 
fireball  and  an  ionized  fireball. 


SPLIT2 

STEVNT 

STOCLD 

STOSET 

TATR 

TILT 

TIMERR 

TIMEW 

TIMSTP 

TRITRP 

TSTCLD 

TVARF2 

TVEF1 

UFUNC 

UPDAT1 

UPDAT2 


Splits  a  mixed  (or  ionized)  fireball  into  a  late  plume 
and  a  subsiding  debris  cloud. 

Routine  to  store  events  in  the  event  list. 

Routine  to  store  cloud  data  in  the  cloud  list. 

Dynamic  storage  routine  to  store  a  D.S.  dataset. 

A  function  used  in  RADC  model  of  height  of  maximum 
electron  density  in  the  F2  layer. 

Used  by  the  propagation  model  to  compute  the  ionospheric 
tilts  due  to  gradients  in  the  ionospheric  electron 
density. 

Logical  function  used  by  the  input  module  to  check  for 
correct  DHMS  format. 

Computes  the  location  in  the  ionospheric  data  table 
and  the  interpolation  weights  for  a  given  local  time. 

Determines  the  appropriate  time  step  for  updating 
mixed,  ionized  or  neutral  fireballs. 

Interpolation  routine  used  by  the  RADC  model  to  obtain  K 
(corrected  Kp)  as  a  function  of  Kp  and  local  geomag¬ 
netic  coordinates;  called  by  SETKP. 

Tests  to  see  if  a  particular  cloud  is  currently  active. 

Function  used  by  the  Aerospace  ionospheric  model. 

Function  used  by  the  Aerospace  ionospheric  model. 

Function  used  by  the  propagation  models. 

UFUNC  =  tan'1  (u)/u-l . 

Handles  the  updating  of  mixed  fireballs  (cloud  type  1). 
Handles  the  updating  of  ionized  fireballs  (cloud  type  2). 
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U PD AT 3 
UPDAT4 
UPDAT5 

VERSY 

VIRTUL 

W 

WIND 

XYMAP 

XYMAPI 

YON  1 1 
ZCGMCS 

ZEUS 

ZFIX 


Handles  the  updating  of  neutral  fireballs  (cloud  type  3) 

Handles  the  updating  of  late  plumes  (cloud  type  4) . 

Handles  the  updating  of  subsiding  debris  clouds  (cloud 
type  5) . 

Evaluates  series  expansions  for  world  maps  of  various 
ITS  ionospheric  parameters. 

Used  by  the  propagation  model  to  compute  virtual  re¬ 
flection  heights. 

Function  used  in  the  Aerospace  ionospheric  model. 

Computes  wind  velocities  for  use  by  the  late-time  sub¬ 
siding  debris  model.  Uses  constant  winds  or  wind 
database  values  averaged  over  the  vertical  thickness 
of  the  cloud  and  interpolated  to  the  latitude  and 
longitude  of  the  cloud  center. 

Converts  latitude  and  longitude  to  map  X-Y  coordinates 
using  the  polyconic,  gnomic  or  "flat  earth"  projection 
equations . 

Converts  map  X-Y  coordinates  to  latitude  and  longitude 
using  the  polyconic,  gnomic  or  "flat  earth"  inverse 
projection  equations. 

Function  used  by  the  Aerospace  ionospheric  model. 

Used  by  the  RADC  (polar)  ionospheric  model  to  convert 
from  geographic  to  (corrected)  geomagnetic  coordinates. 

Computes  atmospheric  noise  power  transmitted  at  a  given 
frequency  from  a  point  source  on  the  equator. 

Adjusts  HF  ray  path  break  points  for  variations  in 
the  virtual  reflection  heights  at  each  hop. 
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MISSION  RESEARCH  CORP  SANTA  BARBARA  CA  '  _  . 

^79*  m°?R0LL1  PR0SR*M  T°  C*LCUL*TE  **'-«*  EFFECTS  On  MF/VMF.i™ 
MRC-N-S1S-V0L-2  DNA-SUTT-1  ‘*«00,-T*-C-0<j2, 


LIST  OF  HFNET  COMMON  BLOCKS 


COMMON 

AMBINT 

BNOISE 

CHEM 

CLOUD 

CON 

CONST 

DATA1 

DEBRIS 


DESCRIPTION  (SIZE) 

Holds  ambient  atmospheric  and  magnetic  properties  at 
the  location  of  a  cloud  before  it  is  updated.  Set  by 
M0D6  and  used  by  all  of  the  cloud  update  routines 
during  a  cloud  update  event.  (11) 

Buffer  for  noise  calculation  dynamic  storage  datasets. 
Used  by  the  propagation  models  in  computing  noise  power 
and  signal -to-noise  ratio.  (50) 

Holds  various  chemistry  constants  and  rate  coefficients. 
Set  by  block  data  routine  BLKDAT  and  used  by  the  D-region 
absorption  models.  (36) 

Holds  cloud  list  bookkeeping  parameters  and  the  cloud 
list  dynamic  storage  dataset  buffer.  Partially  set 
by  block  data  routine  BLKDAT.  (59) 

Holds  various  constants.  Set  by  subroutine  SETPAR 
and  used  by  the  ITS  ionospheric  model  subroutines.  (9) 

Holds  various  constants.  Set  by  block  data  routine 
BLKDAT  and  used  by  many  subroutines.  (8) 

Holds  world  map  data,  interpolated  for  sunspot  number, 
used  to  generate  F2  layer  parameters  in  RADC  ionosphere 
model.  (1019) 

Holds  various  debris  model  parameters  used  in  the  cre¬ 
ation  and  updating  of  clouds.  Set  by  block  data  routine 
BLKDAT.  (31) 

Holds  data  used  by  the  RADC  (polar)  ionospheric  model. 

(9) 


ELAYER 


EVENT 


F2DC0M 

FILES 

GEOPAR 

HDEBUG 

HFLOSS 

HFPATH 

IONBUF 

ITSCOM 

KEYALT 


Holds  event  bookkeeping  parameters  and  3  event  list 
dataset  buffers.  Partially  set  by  block  data  routine 
BLKDAT  and  used  by  the  simulation  manager  and  the 
various  modules.  (52) 

Holds  data  used  to  calculate  upper  and  lower  decile 
frequencies.  The  data  is  read  from  file  DECILE.BIN 
by  subroutine  DPROB  and  is  used  by  subroutine  F2DIS. 
(576) 

Holds  logical  unit  numbers  for  all  input,  output  and 
scratch  files.  Set  by  block  data  routine  BLKDAT  and 
used  throughout  the  code.  (15) 

Holds  various  geographic  and  geomagnetic  parameters 
used  by  the  RADC  (polar)  ionospheric  model .  (9) 

Holds  the  names  of  subroutines  from  which  debug  output 
has  been  req%3Sted.  Initially  set  by  block  data  rou¬ 
tine  BLKDAT  and  possibly  modified  by  input  (using  the 
$DEBUG  keyword) .  (20) 

Holds  data  used  in  the  propagation  loss  calculation 
for  HF  skywave  or  plume  modes.  (587) 

Holds  data  used  in  the  propagation  geometry  calculation 
for  HF  skywave  modes.  Partially  initialized  by  block 
data  routine  BLKDAT.  (369) 

Used  to  buffer  link  ionospheric  data  tables  for  the 
propagation  models.  Partially  initialized  by  block 
data  routine  BLKDAT.  (1591) 

Holds  data  used  by  the  ITS  ionospheric  model.  (3514) 

Holds  data  used  in  the  integration  of  incremental 
D-region  absorption.  Set  by  block  data  routine  BLKDAT 
and  used  in  various  absorption  calculations  (HF  skywave, 
plume  modes,  noise  absorption).  (13) 
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LINK 


Holds  link  list  dataset  buffer.  Partially  set  in 
block  data  routine  BLKDAT.  (11) 

LISTOR  Holds  dynamic  storage  bookkeeping  data  and  the  main 
dynamic  storage  buffer.  Partially  initialized  by 
subroutine  LISTOR  and  used  by  all  of  the  dynamic 
storage  subroutines.  (50099) 

MAPCOM  Holds  map  projection  information  used  by  the  mapping 

routines .  (7) 

MODTRN  Holds  the  module  transfer  list.  Set  by  block  data 

routine  BLKDAT  and  used  by  the  simulation  manager 
(SIMMGR) .  (203) 

MPATH  Holds  data  used  by  the  multimode  propagation  routines. 

(100) 

NEWCLD  Used  to  temporarily  hold  time  variant  cloud  data  for 

a  cloud  that  is  being  created,  updated  or  split.  (20) 

NEWFIX  Used  to  temporarily  hold  time  invariant  cloud  data  for 

a  cloud  that  is  being  created,  updated  or  split.  (10) 

NODE  Holds  node  list  dataset  buffer.  Partially  set  by 

block  data  routine  BLKDAT.  (41) 

NOWDTA  Holds  interpolated  cloud  data  for  a  particular  cloud 
at  a  particular  time.  Set  by  subroutine  INTERP.  (31) 

OLDCLD  Used  to  temporarily  hold  time  variant  data  for  a  cloud 

that  is  being  updated  or  split.  (20) 

OLDFIX  Used  to  temporarily  hold  time  invariant  data  for  a 
cloud  that  is  being  updated  or  split.  (10) 

PLTGRD  Holds  contour  plot  grid  data.  Set  and  used  by  M0D22 
during  contour  plot  events.  (5104) 
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PLUME 


PLUMOD 


POSTER 


RUNTIM 


SUBSID 


TABLIN 


WEATHR 


WHEN 


WHERE 


Holds  plume  parameters .  Set  by  block  data  routine 
BLKDAT  and  used  by  UPDAT4  during  late  plume  update 
events.  (14) 

Holds  data  used  in  plume  mode  propagation  calculations. 
Partially  set  by  block  data  routine  BLKDAT  and  used 
as  a  buffer  for  plume  list  datasets.  (117) 

Holds  data  used  to  generate  a  file  for  input  to  the 
post  processor  program,  HFPOST.  Partially  set  in 
block  data  routine  BLKDAT.  (49) 

Holds  event  and  module  run-time  statistics.  Initialized 
by  block  data  routine  BLKDAT,  updated  by  the  simulation 
manager,  and  printed  by  subroutine  PSTATS.  (202) 

Holds  constants  used  while  updating  subsiding  debris 
clouds.  Set  by  block  data  routine  BLKDAT  and  used  by 
subroutine  UPDAT5.  (6) 

Used  as  an  input  buffer  for  data  read  from  file  GEOMAG.BIN 
and  also  as  a  dynamic  storage  dataset  buffer  for  the 
F77X  data  list.  This  data  is  used  by  the  RADC  (polar) 
ionospheric  model  (subroutines  PREPAR  and  QFETCH) .  (258) 

Used  to  hold  ambient  "weather"  conditions.  Initialized 
by  block  data  routine  BLKDAT  but  usually  reset  by  input 
(using  the  $AMBIENT  keyword) .  Used  throughout  the 
code.  (10) 

Holds  date  and  time  information  for  RADC  (polar) 
ionospheric  model  routines.  Set  by  subroutine  Q0LI0N. 

(8) 

Holds  locational  information  for  RADC  (polar)  ionospheric 
model  routines.  Set  by  subroutine  Q0LI0N.  (7) 
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APPENDIX  C 
IMPLEMENTATION  NOTES 


HFNET  is  written  in  VAX- 11  FORTRAN  IV  PLUS  and  runs  on  a 
DEC  VAX-11/780  minicomputer  using  the  VAX/VMS  operating  system.  Previous 
versions  of  HFNET  have  been  implemented  on  the  following  computer  systems: 

DEC  PDP  11/45  using  RSX-11D, 

DECSYSTEM  10  using  TOPS-10, 

CDC  6600  using  NOS/BE,  and 
CDC  CYBER/176  using  NOS/BE. 

While  most  of  the  code  is  written  in  ANSI  standard  FORTRAN  IV  some  conversion 
would  be  required  to  implement  HFNET  on  another  computer  system.  Conversion 
to  another  VAX,  a  DECSYSTEM  10  or  20,  or  a  large  CDC  system  should  be  a 
straightforward  task,  conversion  to  a  computer  system  significantly  different 
from  these  could  prove  to  be  a  non-trivial  task. 

Some  of  the  known  pitfalls  to  watch  out  for  while  implementing 

HFNET  are: 

•  VAX  machine-dependent  routines:  CPUSEC,  REMTABS,  0PENFILES, 
etc.  These  routines  can  either  be  removed  or  replaced  with 
equivalent  routines. 

•  FORTRAN  77  extensions: 

•  List-directed  input  format 

•  "ERR="  and  "END»"  error  condition  branching 

•  CHARACTER  variables  and  constants 

•  CHARACTER  constants  in  FORMAT  statements 


For  the  most  part  the  uses  of  these  extensions  are  either 
isolated  or  easily  "fixed"  using  a  good  text  editor. 

Some  routines  make  the  assumption  that  arguments  to  sub¬ 
routines  are  passed  "by  reference"  (and  not  "by  value"  or 
some  other  way).  This  is  particularly  true  of  calls  to  the 
array  manipulation  routines  (AMOVX,  ASETX  and  ACLRX) . 

INTEGER  and  REAL  variables  must  occupy  the  same  amount  of 
storage.  This  assumption  is  made  throughout  the  code  in 
conjunction  with  EQUIVALENCE  statements. 

Both  INTEGER  and  REAL  variables  are  used  to  hold  Hollerith 
data.  The  code  assumes  that  4  characters  can  fit  in  a  word 
and  8  into  a  double  precision  word. 

The  dynamic  storage  (D.S.)  buffer  is  dimensioned  to  50,000 
words  with  no  secondary  storage  backup.  The  D.S.  routine 
MEMLOC  can  easily  be  altered  to  allow  for  secondary  storage 
(disk,  large  core  or  whatever  is  available).  This  allows 
the  in-core  D.S.  buffer  to  be  reduced  in  size  without 
restricting  the  total  D.S.  space  available. 

HFNET  is  not  small  -  it  uses  over  500,000  bytes  (including 
the  200KB  D.S.  buffer)  on  the  VAX.  If  this  is  too  large  then 
here  are  3  suggestions: 

1.  Change  the  D.S.  buffersize  (see  above), 

2.  Delete  unwanted  modules  or  models,  or 

3.  Use  an  overlay  scheme. 

HFNET  and  HFPOST  together  require  19  external  files  as  shown 
in  Table  C.l.  How  files  are  implemented  is  very  operating 
system  dependent.  In  order  to  get  an  idea  of  how  they  might 
be  handled,  the  VAX  implementation's  command  procedure  for 
running  the  code  is  shown  in  Table  C.2. 
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Table  C.l.  HFNET/HFPOST  external  files. 


Logical 

Name* 

Internal 

Name 

External 

Name* 

File 

Type** 

File  Description 

**HFNET** 

F0R001 

LUNOAT 

*.DAT 

C/I 

Simulation  input 

F0R002 

LUNOUT 

*.0UT 

P/0 

Detailed  output 

F0R002 

LUNBUG 

*.0UT 

P/0 

Debug  output 

F0R010 

LUNSCN 

*.SCN 

C/I 

Scenario  Input 

F0R012 

LUNF77 

GEOMAG.BIN 

B/I 

RADC  geomagnetic  data 

F0R014 

LUNYR 

YEAR.BIN 

B/I 

ITS  yearly  ionospheric  data 

F0R016 

LUNMON 

M0NTH.BIN 

B/I 

ITS  monthly  ionospheric  data 

F0R020 

LUN20 

*.020 

B/0 

Cloud  plot  data 

F0R021 

LUN21 

*.021 

B/0 

Cloud  list  binary  output 

F0R022 

LUN22 

*.022 

B/0 

Contour  plot  data 

F0R023 

LUN23 

*.023 

B/I 

Cloud  list  binary  input 

F0R050 

LUNPST 

*.PST 

B/0 

Post-processor  data 

F0R060 

LUN1 

*.060 

B/S 

Input  module  scratch  file 

F0R070 

LUNPRB 

DECILE.BIN 

B/I 

ITS  decile  frequency  data 

F0R080 

LUNUND 

WIN0.BIN 

B/I 

Wind  data  base 

**HFP0ST** 

F0R010 

LUNOUT 

SYS$OUTPUT 

P/0 

HFPOST  error  messages 

F0R010 

LUNBUG 

SYS$OUTPUT 

P/0 

HFPOST  debug  output 

F0R020 

LUNPST 

*.  PST 

B/I 

Post-processor  data 

F0R030 

LUNSUM 

*.  SUM 

P/0 

Summary  output 

*  Logical  and  external  file  names  are  for  the  VAX  implementation  of  HFNET. 
**  File  types:  C  =  card  image,  P  *  print,  B  =  binary 
I  =  input,  0  =  output,  S  =  scratch. 
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Table  C.2.  Command  procedure  for  running  HFNET. 


$ 

%  !  HFNETX.COM 
$ 

S  !  VAX/VMS  COMMAND  PROCEDURE  TO  RUN  HFNET  AND  HFP05T 
$ 

$  !  THE  JOB  NAME  IS  PASSED  TO  THIS  PROCEDURE  AS  PARAMETER  1  (PI) 
$ 

$  SET  WORK ING^SET  /LIMIT=200 
* 

■$  JOB^NAME  :=  'PI' 

•i 

$  ASSIGN  ' JOB^NAME' . DAT  FOR001 
$  ASSIGN  ' JOB^NAME' . OUT  FOR002 
$  ASSIGN  ' JOB^NAME' . SCN  FOR010 
$  ASSIGN  6E0MAG.BIN  FOR012 

*  ASSIGN  YEAR .BIN  FOR014 

*  ASSIGN  MONTH .BIN  FOR016 

*  ASSIGN  ’ JOBfNAME' .020  FOR020 
$  ASSIGN  ’ JOB^NAME' .021  FOR021 
%  ASSIGN  ’ JOB*NAME' .022  FOR022 
%  ASSIGN  ' JOB*NAME' .023  FOR023 
%  ASSIGN  ’ JOB^NAME' . PST  FOR050 

*  ASSIGN  ’JOB* NAME' .060  FOR0S0 

$  ASSIGN  DECILE.BIN  FOR070 

$  ASSIGN  UIND.BIN  FOR080 

% 

%  PUN  HFNET 
$ 

$  ASSIGN  SYS$OUTPUT  FOR010  !  TROUBLE  /  DEBUG  OUTPUT 

$  ASSIGN  ’ JOB*NAME' .PST  FOR020  !  POSTPROCESSOR  DATA  FILE 
t  ASSIGN  'JGB^NAME' .SUM  FOR030  !  SUMMARY  OUTPUT 
$ 

$  RUN  HFPOST 
% 

%  DELETE  ' JOB^NAME' .060;* 

%  DELETE  '  JOBc-NAME'  .PST;* 

$ 

$  PRINT  'JOBfNAME' .OUT, 'JOB* NAME' .SUM 

* 

$  SHOU  PROCESS  /ACCOUNTING 
t 


SIMULATION  INPUT 
DETAILED  OUTPUT 
SCENARIO  INPUT 
RADC  GEOMAGNETIC  DATA 
ITS  YEARLY  DATA 
ITS  MONTHLY  DATA 
CLOUD  PLOT  OUTPUT 
CLOUD  LIST  OUTPUT 
CONTOUR  PLOT  OUTPUT 
CLOUD  LIST  INPUT 
POST  PROCESSOR  DATA  FILE 
INPUT  MODULE  SCRATCH  FILE 
ITS  DECILE  DATA 
UIND  DATA 


The  code  comes  on  a  standard  9-track,  unlabeled,  800  BPI,  ASCII 
magnetic  tape.  This  tape  contains  the  following  11  files: 

1 .  HFNET  program 

2.  HFPOST  program 

3.  A  program  (CONVERT)  to  convert  the  ASCII  data  files 
(files  4-7)  to  binary  format  for  use  by  HFNET 

4.  ITS  monthly  data  file 

5.  ITS  yearly  data  file 

6.  RADC  geomagnetic  data  file 

7.  Wind  data  file 

8.  VAX/ VMS  command  procedure  for  setting  up  the  program 

9.  VAX/ VMS  command  procedure  for  running  the  program 

10.  Test  case  data  file 

11.  Test  case  output  file 

The  files  on  the  tape  are  all  in  FIXED  RECORD  SIZE/BLOCKED  format  with  the 
record  and  block  sizes  (in  bytes)  as  shown  in  Table  C.3  below. 

Table  C.3.  HFNET  tape  file  attributes. 


File 

# 

File  Name 

Record 

Size 

Block 

Size 

# 

Records 

# 

Blocks 

File 

Type 

1 

HFNET. FOR 

88 

4400 

26583 

532 

source 

2 

HFPOST. FOR 

88 

4400 

2443 

49 

II 

3 

CONVERT. FOR 

80 

4000 

138 

3 

II 

4 

M0NTH.DAT 

80 

4000 

24660 

494 

data 

5 

YEAR.DAT 

120 

3600 

4006 

134 

II 

6 

GE0MAG.DAT 

80 

4000 

1314 

27 

II 

7 

WIND.DAT 

80 

4000 

7884 

158 

II 

8 

SETUP.COM 

80 

4000 

18 

1 

command 

9 

HFNETX.COM 

80 

4000 

41 

1 

II 

10 

TEST.DAT 

80 

4000 

72 

2 

data 

11 

TEST. LIS 

133 

2660 

1056 

53 

print 

DISTRIBUTION  LIST 
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Defense  Communications  Engineer  Center 

ATTN 

Code  R410,  N.  Jones 

ATTN:  Code  R123 

Defense  Intelligence  Agency 

ATTN 

DB-4C ,  E.  O' Farrell 

ATTN 

HQ-TR,  J.  Stewart 

ATTN 

DB,  A.  Wise 

ATTN 

DT-5 

ATTN 

DT-1B 

ATTN:  DC-70,  W.  Wittig 

Defense  Nuclear  Agency 

ATTN 

RAEE 

ATTN 

STNA 

ATTN 

NATD 

ATTN 

NAFD 

4  cy  ATTN 

TITL 

3  cy  ATTN 

RAAE 

Defense  Technical  Information  Center 

12  cy  ATTN 

DD 

Field  Command 
Defense  Nuclear  Agency 
ATTN:  FCPR 


Field  Command 

Defense  Nuclear  Agency 
ATTN:  FCPRL 

Interservice  Nuclear  Weapons  School 
ATTN:  TTV 

Joint  Chiefs  of  Staff 

ATTN:  C3S  Evaluation  Office 
ATTN:  C3S 

Joint  Strat  Tgt  Planning  Staff 
ATTN:  JLTW-2 
ATTN:  JLA 

Under  Secy  of  Def  for  Rsch  I  Engrg 

ATTN:  Strategic  4  Space  Sys  (OS) 


DEPARTMENT  OF  DEFENSE  (Continued) 

National  Security  Agency 

ATTN:  W-32,  0.  Bartlett 
ATTN:  B-3,  F.  Leonard 
ATTN:  R-52,  J.  Skillman 

WWMCCS  System  Engineering  Org 
ATTN:  J.  Hoff 

DEPARTMENT  OF  THE  ARMY 

Assistant  Chief  of  Staff  for  Automation  4  Conrn 
Department  of  the  Army 

ATTN:  DAAC-ZT,  P.  Kenny 

Atmospheric  Sciences  Laboratory 
U.S.  Army  Electronics  R4D  Comand 
ATTN:  DELAS-EO,  F.  Niles 

BMD  Advanced  Technology  Center 
Department  of  the  Army 

ATTN:  ATC-T,  M.  Capps 
ATTN:  ATC-0,  W.  Davies 

BMD  Systems  Command 
Department  of  the  Army 
2  cy  ATTN:  BMOSC-HW 

Dep  Ch  of  Staff  for  Ops  4  Plans 
Department  of  the  Army 
ATTN:  DAMO-RQC 

Electronics  Tech  4  Devices  Lab 
U.S.  Army  Electronics  R4D  Command 
ATTN:  DELET-ER,  H.  Bomke 

Harry  Diamond  Laboratories 
Department  of  the  Army 

ATTN:  DELHD-N-P,  F.  Wimenitz 
ATTN:  DELHD-N-P 
ATTN:  DELHD-N-RB,  R.  Williams 
ATTN:  DELHD-I-TL,  M.  Weiner 

U.S.  Army  Conm-Elec  Engrg  Instal  Agency 
ATTN:  CCC-CED-CCO,  W.  Neuendorf 
ATTN:  CCC-EMEO-PED,  G.  Lane 

U.S.  Army  Communications  Command 
ATTN:  CC-OPS-W 
ATTN:  CC-OPS-WR,  H.  Wilson 

U.S.  Army  Conimmlcatlons  R4D  Command 

ATTN:  DRDCO-COM-RY,  W.  Kesselman 

U.S.  Army  Foreign  Science  4  Tech  Ctr 
ATTN:  DRXST-SD 

U.S.  Army  Materiel  Dev  4  Readiness  Ond 
ATTN:  DRCLDC,  J.  Bender 

U.S.  Army  Missile  Intelligence  Agency 
ATTN:  J.  Gamble 

U.S.  Army  Nuclear  4  Chemical  Agency 
ATTN:  Library 


DEPARTMENT  OF  THE  ARMY  (Continued) 

U.S.  Army  Satellite  Comm  Agency 
ATTN:  Document  Control 

U.S.  Army  TRAOOC  Sys  Analysis  Actvy 
ATTN:  ATAA-TCC,  F.  Payan,  Jr 
ATTN:  ATAA-TDC 
ATTN:  ATAA-PL 

DEPARTMENT  OF  THE  NAVY 

Joint  Cruise  Missiles  Project 

Department  of  the  Navy 
ATTN:  JCMG-707 

Naval  Air  Development  Center 

ATTN:  Code  6091,  M.  Setz 

Naval  Air  Systems  Command 
ATTN:  PMA  271 

Naval  Electronic  Systems  Command 

ATTN:  PME-117-2013,  G.  Burnhart 

ATTN :  Code  501A 

ATTN:  Code  3101,  T.  Hughes 

ATTN:  PME  106-13,  T.  Griffin 

ATTN:  PME  117-20 

ATTN:  PME  106-4,  S.  Kearney 

ATTN:  PME  117-211,  B.  Kruger 

Naval  Intelligence  Support  Ctr 
ATTN:  NISC-50 

Naval  Ocean  Systems  Center 


ATTN: 

Code 

532, 

J. 

Bickel 

ATTN: 

Code 

5322, 

M. 

Paulson 

3  cy 

ATTN: 

Code 

5324, 

W. 

Moler 

3  cy 

ATTN: 

Code 

5323, 

J. 

Ferguson 

Naval 

Research  Laboratory 

ATTN: 

Code 

4780, 

S. 

Ossakow 

ATTN: 

Code 

7550, 

J. 

Davis 

ATTN: 

Code 

4187 

ATTN: 

Code 

4700, 

T. 

Coffey 

ATTN: 

Code 

7500, 

B. 

Wald 

Naval  Space  Surveillance  System 
ATTN:  J.  Burton 


Naval  Surface  Weapons  Center 
ATTN:  Code  F31 

Naval  Telecommunications  Command 
ATTN:  Code  341 

Office  of  Naval  Research 
ATTN:  Code  420 
ATTN:  Code  465 
ATTN:  Code  421 

Office  of  the  Chief  of  Naval  Operations 
ATTN:  OP  65 
ATTN:  OP  98 IN 
ATTN:  OP  9410 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
ATTN:  NSP-2141 
ATTN:  NSP-43 

ATTN:  NSP-2722,  F.  Wimberly 


DEPARTMENT  OF  THE  AIR  FORCE 

Aerospace  Defense  Command 
Department  of  the  Air  Force 
ATTN:  DC,  T.  Long 

Air  Force  Geophysics  Laboratory 
ATTN:  OPR,  H.  Gardiner 
ATTN:  OPR-1 
ATTN:  LKB,  K.  Champion 
ATTN:  OPR,  A.  Stair 
ATTN:  PHP 

ATTN:  PHI,  J.  Buchau 

Air  Force  Weapons  Laboratory 
Air  Force  Systems  Command 
ATTN:  SUL 
ATTN:  NTYC 
ATTN:  NTN 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  W.  Hunt 
ATTN:  A.  Johnson 

Air  Logistics  Conmand 
Department  of  the  Air  Force 

ATTN:  OO-ALC/MM,  R.  Blackburn 

Air  University  Library 
Department  of  the  Air  Force 
ATTN:  AUL-LSE 

Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
ATTN:  INED 

Assistant  Chief  of  Staff 
Studies  A  Analyses 
Department  of  the  Air  Force 

ATTN:  AF/SASC,  G.  Zank 
ATTN:  AF/SASC,  W.  Adams 

Ballistic  Missile  Office 
Air  Force  Systems  Command 

ATTN:  MNNXH ,  J.  Allen 

Deputy  Chief  of  Staff 
Operations  Plans  and  Readiness 
Department  of  the  Air  Force 
ATTN:  AFXOKT 
ATTN:  AFXOXFD 
ATTN:  AFXOKCD 
ATTN:  AFXOKS 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq 
Department  of  the  Air  Force 
ATTN:  AFRDSP 
ATTN:  AFRDS 
ATTN:  AFRDSS 
ATTN:  AFRDS 

Electronic  Systems  Division 
ATTN:  DCKC,  J.  Clark 

Electronic  Systems  Division 
ATTN:  XRW,  J.  Deas 
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DEPARTMENT  Of  THE  AIR  FORCE  (Continued) 

Electronic  Systems  Division 
ATTN:  YSM,  J.  Kobelski 
ATTN:  YSEA 

Foreign  Technology  Division 
Air  Force  Systems  Command 

ATTN:  TQTD,  B.  Ballard 
ATTN:  NIIS  Library 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  SKA,  D.  Bolin 
ATTN:  SKY,  C.  Kennedy 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  YZJ ,  W.  Mercer 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN,-  E.  Butt 

Rome  Air  Development  Center 
Air  Force  Systems  Conmand 
ATTN:  OCS,  V.  Coyne 
ATTN:  TSLD 

Rome  Air  Development  Center 
Air  Force  Systems  Command 
ATTN:  EEP 

Strategic  Air  Conmand 
Department  of  the  Air  Force 

ATTN:  DCXR,  T.  Jorgensen 

ATTN:  DC  XT 

ATTN:  NRT 

ATTN:  XPFS 

ATTN:  DCX 

OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NEO 

Department  of  Commerce 
National  Bureau  of  Standards 

ATTN:  Sec  Ofc  for  R.  Moore 

Department  of  Commerce 
National  Oceanic  &  Atmospheric  Admin 
ATTN:  R.  Grubb 

Institute  for  Telecommunications  Sciences 
ATTN:  H.  Utlaut 
ATTN:  A.  Jean 
ATTN:  L.  Berry 

U.S.  Coast  Guard 
Department  of  Transportation 

ATTN:  G-DOE-3/TP54,  B.  Romine 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Lawrence  Livermore  National  Lab 
ATTN:  L-31,  R.  Hager 
ATTN:  Technical  Info  Dept  Library 
ATTN:  L-389,  R.  Ott 


DEPARTMENT  OF  ENERGY  CONTRACTORS  (Continued) 

Los  Alamos  National  Scientific  Lab 
ATTN:  MS  670,  J.  Hopkins 
ATTN:  MS  664,  J.  Zinn 
ATTN:  D.  Simons 
ATTN:  E.  Jones 
ATTN:  D.  Westervelt 
ATTN:  P.  Keaton 
ATTN:  R.  Taschek 

Sandia  Laboratories 

Livermore  Laboratory 
ATTN:  T.  Cook 
ATTN:  B.  Murphey 

Sandia  National  Lab 

ATTN:  Org  1250,  W.  Brown 
ATTN:  D.  Thornbrough 
ATTN:  D.  Dahlgren 
ATTN:  Space  Project  Div 
ATTN:  3141 

ATTN:  Org  4241,  T.  Wright 

EG&G,  Inc 

ATTN:  J.  Colvin 
ATTN:  D.  Wright 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  I.  Garfunkel 
ATTN:  D.  Olsen 
ATTN:  N.  Stockwell 
ATTN:  R.  Slaughter 
ATTN:  S.  Bower 
ATTN:  J.  Straus 
ATTN:  V.  Josephson 
ATTN:  T.  Salmi 

University  of  Alaska 
ATTN:  N.  Brown 
ATTN:  T.  Davis 
ATTN:  Technical  Library 

Analytical  Systems  Engineering  Corp 
ATTN:  Radio  Sciences 

Analytical  Systems  Engineering  Corp 
ATTN:  Security 

Barry  Research  Corporation 
ATTN:  J.  McLaughlin 

BDM  Corp 

ATTN:  T.  Neighbors 
ATTN:  L.  Jacobs 

Berkeley  Research  Associates,  Inc 
ATTN:  J.  Workman 

BETAC 

ATTN:  J.  Hirsch 

Boeing  Co 

ATTN:  G.  Hall 

ATTN:  S.  Tashird 

ATTN:  M/S  42-33,  J.  Kennedy 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Booz -Allen  &  Hamilton,  Inc 

Johns  Hopkins  University 

ATTN:  B.  Wilkinson 

ATTN:  P.  Komiske 

University  of  California  at  San  Diego 

ATTN:  T.  Potemra 

ATTN:  T.  Evans 

ATTN:  H.  Booker 

ATTN:  J.  Newland 

ATTN:  J.  Phillips 

Charles  Stark  Draper  Lab,  Inc 

ATTN:  D.  Cox 

Kaman  TEMPO 

ATTN:  0.  Gilmore 

ATTN:  W.  Knapp 

Computer  Sciences  Corp 

ATTN:  M.  Stanton 

ATTN:  T.  Stephens 

ATTN:  H.  Blank 

ATTN:  DAS  I AC 

ATTN:  W.  McNamara 

Comsat  Labs 

ATTN:  G.  Hyde 

ATTN:  R.  Taur 

Linkabit  Corp 

ATTN:  I.  Jacobs 

Cornell  University 

ATTN:  D.  Farley,  Jr 

Litton  Systems,  Inc 

ATTN:  M.  Kelly 

ATTN:  R.  Grasty 

Electrospace  Systems,  Inc 

Lockheed  Missiles  &  Space  Co, 

Inc 

ATTN:  H.  Logston 

ATTN:  W.  Imhof 

ESL,  Inc 

ATTN:  R.  Johnson 

ATTN:  M.  Walt 

ATTN:  J.  Marshall 

Lockheed  Missiles  &  Space  Co, 

Inc 

General  Electric  Co 

ATTN:  D.  Churchill 

ATTN:  A.  Ha  rear 

ATTN:  Dept  60-12 

ATTN:  M.  Bortner 

General  Electric  Co 

M.I.T.  Lincoln  Lab 

ATTN:  L.  Loughlin 

ATTN:  C.  Zierdt 

ATTN:  D.  Towle 

ATTN:  A.  Steinmayer 

General  Electric  Co 

Martin  Marietta  Corp 

ATTN:  R.  Heffner 

ATTN:  F.  Reibert 

General  Electric  Tech  Services  Co,  Inc 

McDonnell  Douglas  Corp 

ATTN:  J.  Moule 

ATTN:  G.  Millman 

ATTN:  W.  Olson 

General  Research  Corp 

ATTN:  G.  Mroz 

ATTN:  R.  Halprin 

ATTN:  J.  Garbarino 

ATTN:  N.  Harris 

ATTN:  J.  Ise,  Jr 

Meteor  Communications  Consultants 

HSS,  Inc 

ATTN:  R.  Leader 

ATTN:  D.  Hansen 

IBM  Corp 

Mission  Research  Corp 

ATTN:  R.  Bogusch 

ATTN:  F.  Ricci 

ATTN:  Tech  Library 

University  of  Illinois 

ATTN:  D.  Sappenfield 
ATTN:  R.  Kilb 

ATTN:  K.  Yeh 

ATTN:  S.  Gutsche 

Institute  for  Defense  Analyses 

ATTN:  F.  Fajen 

ATTN:  R.  Hendrick 

ATTN:  J.  Bengston 

5  cy  ATTN:  Dociment  Control 

ATTN:  H.  Wolfhard 

ATTN:  J.  Aein 

Mitre  Corp 

ATTN:  E.  Bauer 

ATTN:  B.  Adams 

International  Tel  &  Telegraph  Corp 

ATTN:  C.  Callahan 

ATTN:  G.  Harding 

ATTN:  G.  Wetmore 

ATTN:  A.  Kytimel 

ATTN:  Technical  Library 

JAYCOR 

Mitre  Corp 

ATTN:  M,  Horrocks 

ATTN:  J.  Sperling 

ATTN:  J.  Wheeler 

JAYCOR 

ATTN:  W.  Foster 

ATTN:  W.  Hall 

ATTN:  J.  Donearlos 
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Pacific-Sierra  Research  Corp 
ATTN:  E.  Field,  Jr 
ATTN:  F.  Thomas 
ATTN:  H.  Brode 

Pennsylvania  State  University 

ATTN:  Ionospheric  Research  Lab 

Photometries,  Inc 


ATTN:  I. 

Kofsky 

Physical  Dynamics,  Inc 

ATTN 

E. 

Fremouw 

R  &  D  Associates 

ATTN 

C. 

Greif inger 

ATTN 

B. 

Gabbard 

ATTN 

R. 

Lelevier 

ATTN 

F. 

Gilmore 

ATTN 

U. 

Karzas 

ATTN 

R. 

Turco 

ATTN 

H. 

Ory 

ATTN 

M. 

Gantsweg 

ATTN 

W. 

Wright 

ATTN 

P. 

Haas 

R  &  D  Associates 

ATTN 

L. 

Delaney 

ATTN 

B. 

Yoon 

Rand  Corp 

ATTN 

E. 

Bedrozian 

ATTN 

C. 

Crain 

Riverside  Research  Institute 

ATTN 

V. 

Trapani 

Rockwell  International  Corp 
ATTN:  R.  Buckner 


Rockwell  International  Corp 
ATTN:  S.  Quilici 

Santa  Fe  Corp 

ATTN:  D.  Paolucci 

Science  Applications,  Inc 


ATTN: 

E. 

Straker 

ATTN: 

D. 

Hamlin 

ATTN: 

L. 

Llnson 

ATTN: 

C. 

Smith 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Science  Applications,  Inc 
ATTN:  SZ 


Science  Applications,  Inc 
ATTN:  J.  Cockayne 


SRI  International 

ATTN: 

R. 

Livingston 

ATTN: 

G. 

Price 

ATTN: 

A. 

Burns 

ATTN: 

C. 

Rino 

ATTN: 

J. 

Petri ekes 

ATTN: 

R. 

Tsunoda 

ATTN: 

M. 

Baron 

ATTN: 

R. 

Leadabrand 

ATTN: 

D. 

Neil  son 

ATTN: 

G. 

Smith 

ATTN: 

W. 

Jaye 

ATTN: 

W. 

Chesnut 

Syl vania  Systems  Group 
ATTN:  M.  Cross 


Technology  International  Corp 

ATTN:  Doc  Con  For  W.  Boquist 

Tri-Com,  Inc 

ATTN:  D.  Murray 

TRW  Defense  &  Space  Sys  Group 
ATTN:  D.  Dee 
ATTN:  R.  Plebuch 

Utah  State  University 
ATTN:  K.  Baker 
ATTN:  L.  Jensen 
ATTN:  J.  Dupnik 

Vi  si dyne,  Inc 

ATTN:  J.  Carpenter 
ATTN:  C.  Humphrey 
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